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Abstract. In this article we study in detail a family of random matrix ensembles which 
are obtained from random permutations matrices (chosen at random according to the Ewens 
f^ ' measure of parameter > 0) by replacing the entries equal to one by more general non- 

04 . vanishing complex random variables. For these ensembles, in contrast with more classical 

j;».^' models as the Gaussian Unitary Ensemble, or the Circular Unitary Ensemble, the eigen- 

C^ . values can be very explicitly computed by using the cycle structure of the permutations. 

Moreover, by using the so-called virtual permutations, first introduced by Kerov, Olshanski 
and Vershik, and studied with a probabilistic point of view by Tsilevich, we are able to de- 
^^ ' fine, on the same probability space, a model for each dimension greater than or equal to one, 

which gives a meaning to the notion of almost sure convergence when the dimension tends 
r^ . ' to infinity. In the present paper, depending on the precise model which is considered, we 

p I . obtain a number of different results of convergence for the point measure of the eigenvalues, 

_5 ' some of these results giving a strong convergence, which is not common in random matrix 

-*— > . theory. 
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1. Introduction 

1.1. Random permutation matrices and outline of the paper. The distribution of 
the eigenvalues of random matrices and some related objects such as their characteristic 
polynomials have received much attention in the last few decades. They have been applied in 
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such diverse branches as physics, number theory, analysis or probabihty theory as illustrated 
in the monographs [11], [10] and [1] or the survey paper f5]. The main matrix ensembles 
which have been studied are the Gaussian ensembles and some of their generalizations, the 
classical compact Lie groups U{N) (the group of unitary matrices), 0{N) (the orthogonal 
group) and some of their continuous sub-groups, endowed with the Haar probability measure. 
It is also natural to investigate for the distribution of the eigenvalues of random permutation 
matrices, i.e. matrices which are canonically associated to a random element of a given finite 
symmetric group. Indeed it is well-known that the eigenvalues of a permutation matrix Mg- 
associated with a permutation a are entirely determined by the cycle structure of a, and 
hence one can hope to take advantage of the extensive literature on random permutations 
(see e.g. the book by Arratia, Barbour and Tavare [2J) to describe completely the structure 
of the point process of the eigenvalues of random permutation matrices (e.g. the correlation 
measure of order q, the convergence of the normalized and non-normalized empirical spectral 
distribution, etc.). This has been shortly sketched out in the pioneering work by Diaconis 
and Shahshahani [6] and further developed by Wieand in [IS], who studied the problem of 
counting how many eigenvalues lie in some fixed arc of the unit circle. Wieand compares 
the results obtained in this case with those obtained for the unitary group under the Haar 
measure and notices some similarities but also some differences when it comes to look at 
more refined aspects. Then it is suggested that one should try to compute finer statistics 
related to the eigenvalues in order to see how random permutation matrices fit into the 
random matrix picture. Of course one expects some drastic differences: for instance the 
point process associated with the eigenvalues of random permutation matrices should not 
be determinantal whereas it is determinantal for the unitary group. The goal of the present 
work is to continue the work initiated by Diaconis and Shahshahani in [B] and Wieand in 
^16j. Before mentioning more precisely the various directions in which we wish to extend the 
existing work, it should be mentioned that other works have been recently done on random 
permutation matrices, such as the paper by Hambly, Keevach, O'Connell and Stark [8] on the 
characteristic polynomial of random unitary matrices, the papers by Zeindler ([H]) and by 
Dehaye and Zeindler ([^), or the works by Wieand ([H]) and Evans ([7J) on the eigenvalues 
of random wreath products (but in this latter case, the techniques that are involved are 
different and we shall also address this framework in a future work). 

We now briefly mention the way we shall continue and extend some of the previous works 
and shall postpone precise definitions to the next paragraph: 



We shall consider a larger ensemble of random matrices with more general distribu- 
tions; roughly speaking, we first pick a permutation of size N from S^r (the group 
of permutations of size N) at random according to the Ewens measure of parameter 
^ > (under this measure, the probability of a given permutation a is proportional 
to 6'", where n is the number of cycles of a) and then consider the corresponding 
matrix; we then replace the I's by a sequence of i.i.d. random variables zi, . . . ,zn 
taking values in C*. This ensemble is a sub-group of the linear sub-group GL{N, C) 
and the classical random permutation matrices correspond to the case where 6 = 1 
and where the distribution of the ^j's is the Dirac measure at 1. The choice of the 
Ewens measure is natural since it is a one parameter deformation of the uniform 
distribution which is coherent with the projections from S^' onto S^v for N > N 
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(see the next paragraphs). If the Zi^s take their values on the unit circle, then our 
ensemble is a sub-group of U{N). 

• We shall also give a meaning to almost sure convergence for the probability empirical 
spectral distribution. Indeed, it is not a priori obvious to define almost sure conver- 
gence in random matrix theory since the probability space changes with N and to 
the best of our knowledge there has not been, so far, a satisfactory way to overcome 
this difficulty. We propose to deal with this problem by considering the space of vir- 
tual permutations, which were introduced by Kerov, Olshanski and Vershik in [9] . A 
virtual permutation is a sequence of permutations (crAr)Ar>i constructed in a coherent 
way, in the sense that for all A^' > A^ > 1, the cycle structure of (Jn can be obtained 
from the cycle structure of a^y' simply by removing the elements strictly larger than 
A^. The virtual permutations satisfy the following property: if (o"Ar)7v>i is a virtual 
permutation, if A^' > A^ > 1 and if a^y' is distributed according to the Ewens mea- 
sure of parameter 6 on S^/ , then a^ follows the Ewens measure of parameter 6 on 
S^r. We shall then take advantage of the work by Tsilevich [H], who proved almost 
sure convergence for the relative cycle lengths of virtual permutations. We shall also 
extensively use an algorithmic way to generate the Ewens measure on the space of 
virtual permutations. 

• We shall study in detail the point process of the eigenvalues of a matrix M drawn 
from our matrix ensembles. For instance we establish various convergence results for 
the empirical spectral measure n{M) = ^ S\, or for fi{M)/N {N being the dimension 
of the matrix) where the sum is over all eigenvalues A of M, counted with multiplicity. 
We also consider the average measure obtained by taking the expectation of fi{M), 
as well as the correlation measure of order q. In the special case where the random 
variables Zi^s take their values on the unit circle, the point process associated with 
the eigenvalues has some more remarkable properties. For instance, when the Zi^s 
are uniformly distributed on the unit circle, the empirical spectral measure as well 
as its limit (in a sense to be made precise) are invariant by translation. Still in this 
case, the 1-correlation for fixed A^ and A^ = oo is the Lebesgue measure, whereas 
for g > 2, and fixed A^, the g-correlation is not absolutely continuous with respect 
to the Lebesgue measure anymore. For A^ = oo, the pair correlation measure is 
still absolutely continuous with respect to the Lebesgue measure with an explicit 
density, but this result fails to hold for the correlations of order greater than or 
equal to 3. One can push further the analogy with the studies made for the classical 
compact continuous groups by characterizing the distribution function of the smallest 
eigenangle (again when the Zi^s are uniformly distributed on the unit circle) as the 
solution of some integral equation. 

1.2. Definitions and notation. In this section we describe the structure of the eigenvalues 
of the "generalized permutation matrices" mentioned above, independently of any probability 
measure, and then we define the family of probability measures which will be studied in detail 
in the sequel of this article. More precisely let Sat be the group of permutations of order 
A^. By straightforward computations, one can check that the set of matrices M such that 
there exists a permutation a G S/y and complex numbers Zi,Z2, ■■■,Zn, different from zero, 
satisfying Mjk = Zjlj=^(^k), is a multiplicative group, denoted by ^(A^) in this article, and 
which can be written as the wreath product of C* and S^y. The group Q^ can also be 
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viewed as the group generated by the permutation matrices and the diagonal matrices of 
GL{N,C). The elements of Q{N) such that \zj\ = 1 for all j, 1 < j < A^, form a subgroup 
noted 'H{N) of Q{N), which can be viewed as the wreath product of U (the set of complex 
numbers with modulus one) and S^v. One can also define, for all integers k > 1, the group 
'Hfc(A^) of elements of ^(A^) such that Zj = I ior I < j < N: the subgroup 'Hfc(A^) of 
'H(A^) is the wreath product of Uk (group of k-th roots of unity) and Sat. Note that the 
group structure oi Q, Ti or T-Ln does not play a fundamental role in our work, however, this 
structure gives a parallel between this paper and the study of other groups of matrices, such 
as the orthogonal group, the unitary group or the symplectic group (we note that "H and Ti^ 
are subgroups of the unitary group of dimension A^). The advantage of the study of matrices 
in ^(A^) is the fact that the structure of their eigenvalues can be very explicitly described in 
function of the cycle structure of the corresponding permutations. More precisely, let M be 
a matrix satisfying Mj^ = -2jlj=o-(fc) for all 1 < j,k < N, where a is a permutation of order 
A^, and Zi, ..., Zat G C*. If the supports of the cycles of a are Ci, ..., C„, with corresponding 
cardinalities /i, ...,/„, and if for 1 < m < n, Rm is the set of the roots of order /„ of the 
complex number 

^m • 11 Zj) 

jeCm 
then the set of eigenvalues of M is the union of the sets Rm, and the multiplicity of any 
eigenvalue is equal to the number of sets Rm containing it. An example of explicit calculation 
is that of the trace: since for all integers / > 2, the sum of the /-th roots of unity is equal to 
zero, one immediately deduces that 

Tr(M)=5^^„ 

where F is the set of fixed points of a. More generally, one can compute the trace of all 
the powers of M. Indeed, for all integers k > 1, the eigenvalues of M^ can be computed by 
taking the k-th powers of the elements of the sets Rm- Therefore, 

n 

Tr(M^) = 5^ Y, ujK 

Now, if Im is not a divisor of k, the last sum is equal to zero, and if /„ is a divisor of k, all 
the terms of the last sum are equal to Zj "". We deduce that 

Tr(M'=) = ^/„Z^/''». 

Lm I f^ 

We see that the description of the eigenvalues and the computations above do not depend 
on any probability measure given on the space Q{N). Let us now define a particular class of 
probability measures on ^(A^) which will be studied in detail in this paper. 

Definition 1.1. Let ^ > and let £ be a probability distribution on C*. The probability 
measure P(A^, 9, C) on ^(A^) is the law of the matrix M{a, Zi, ..., z^), where: 

• the permutation a follows the Ewens measure of parameter 6 on E^r, i.e. the prob- 
ability that 0" is equal to a given permutation is proportional to 6^\ where n is the 
number of cycles of a. 

• for all J, 1 < J < A^, Zj is a random variable following the probability law C. 



• the random permutation a and the random variables zi, ...,zn are all independent. 

• M{a,zi, ...,zn) is the matrix M G ^(A^) such that for all 1 < j,k < N, Mjk = 

In this paper we prove for a large class of probability distributions C, the weak convergence 
of the law of the empirical measure of the eigenvalues, when the dimension N tends to infinity. 
In several particular cases we are also interested in almost sure convergences, and then we 
need to couple all the dimensions A^ on the same probability space. This can be done by 
introducing the so-called virtual permutations, which were first defined by Kerov, Olshanski 
and Vershik in [9] and also studied by Tsilevich [Hj. A virtual permutation is a sequence 
{o'n)n>i of permutations, such that for all N > 1, aN & ^n, and the cycle structure of a^ 
is obtained from the cycle structure of ctat+i, simply by removing the element A^ + 1 (for 
example if ag = (13745)(28)(6), then ^7 = (13745)(2)(6)). Now, for 9 > 0, it is possible 
to define on the space of virtual permutations the so-called Ewens measure of parameter 6 
as the unique (by the monotone class theorem) probability measure, such that if {crN)N>i 
follows this measure, a^ follows the Ewens(6') measure on Stv. Now we can introduce the 
following definition: 

Definition 1.2. Let 6 > and let £ be a probability law on C*. The probability measure 
P(oo,6',£), defined on the product of the probability spaces ^(A^), A^ > 1, is the law of a 
sequence of random matrices (MAr)jv>i, such that M^ = M{aN, zi, ..., z^), where: 

• the sequence {(Tn)n>i is a random virtual permutation following Ewens measure of 
parameter 6. 

• for all J > 1, Zj is a random variable following the distribution C. 

• the virtual permutation (cr7v)Ar>i and the random variables {zj)j>i are independent. 

It is easy to check that for all A^ > 1, the image of F{oo,6,C) by the A^-th coordinate 
projection is the measure P(A^, ^,£). The properties on F{N,9,C) or F{oo,6,C) which are 
obtained in this article depend on the probability distribution C in an essential way. 

In the next section, we shall review some properties of virtual permutations which are 
needed in the sequel of the paper, and then we give general results on the point process 
associated with the eigenvalues of random matrices from Q{N). We finally refine some of 
these results in the case of T-L{N). 

2. Generating the Ewens measure on the set of virtual permutations 

The space of virtual permutations was first introduced by Kerov, Olshanski and Vershik 
in [H] in the context of representation theory; the interested reader can refer to the notes by 
Olshanski in [12] for more details and references. Here we shall mostly be concerned with 
the probabilistic aspects of virtual permutations which were studied in detail by Tsilevich in 
[H] . We now review a construction of the virtual permutations which is explained in [H] and 
which is suitable for the probabilistic reasoning. We then show how to generate the Ewens 
measure of parameter 6 on the space of virtual permutations. This is already explained by 
other means by Tsilevich in her unpublished note [I5]; here we provide a more elementary 
way to generate it and give all details since this is going to be at the source of many of our 
proofs and because also virtual permutations might not be so well-known. In the sequel, we 
shall assume that the reader is familiar with the GEM and Poisson-Dirichlet distributions 
(if not, one can refer to [3], p. 40-48). 
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As already mentioned, one of the interests of the construction of virtual permutations is 
that it gives a natural explanation of convergences in law involved when one looks at the 
relative lengths of cycles of random permutations in the symmetric group S^v when A^ goes 
to infinity. Recall that for all integer A^ > 1, there exists a bijective map 

N 

$^ : J]{l,2,...,j} ^ S;v 
i=i 
such that 

where for Tj^k is the unique permutation such that Tj^k{j) = k, Tj^kik) = j, and Tj^kiO = I 
for / different from j and k (if j ^ k, Tj^k is a transposition, if j = k, it is the identity). If 
A"! > A'2, the bijections ^Ni and ^n2 induce a natural surjective map ttni,N2 from SjVi to 
SjVa, defined in the following way: if a G S^r-^, there exists a unique sequence {mj)i<j<^-^, 
rrij G {1, ..., j} such that 

a = <l>jvi (("^i)i<j<jvj , 
and one then defines: 

Note that for A^i > A'2 > A'3, 

Now, a virtual permutation is a sequence of permutations (criv)jv>i, such that ajy G S^v for 
all A^, and which is consistent with respect to vr, i.e. for all A^i > A"2, 

We denote by Sqo the (infinite) set of virtual permutations, and note that the group structure 
of Sat does not induce any group structure on Eqo. Moreover, there is a natural bijection, 

00 
$00 : l[{l,2,...,j} ^ Soo 
i=i 

induced by the bijections $Ar. Indeed, for all infinite sequences {'mj)j>i, nij G {1, ..., j} one 
defines: 

*oo (("^i)i>i) := (crjv)Af>i 
where for all A^, 

aN = ^N (("^j)i<i<jv) • 
Here, one immediately checks that {aN)N>i is consistent. Moreover, one can also define a 
surjective map iIoo,n from Eoo to E^v for all A^ > 1, by setting: 

TToo.Af ((0'n)n>l) := CTn, 

and one can check the relation 

The cycle structure of a virtual permutation can be described by the so-called Chinese 
restaurant process, described for example by Pitman ([13]). More precisely, let 

O'oo = {o'n)n>1 G Soo 
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be a virtual permutation. There exists a unique sequence ("^j)j>l, rrij G {1, ..., j} such that 

Then the cycle structure for ctat can be described by induction on A^: 

• 0"! is (of course!) the identity of Si; 

• if m7v+i = A^ + 1 for iV > 1, the cycle structure of aN+i is obtained from the structure 
of (Jiv by simply adding the fixed point N + 1; 

• if niN+i < N ioT N > 1, the cycle structure of aN+i is obtained from the structure 
of (Tat by inserting A^ + 1 just before rriN+i, in the cycle containing mTv+i- 

For example, if a^ = (124)(3) and mg = 5, then as = (124)(3)(5), and if a^ = (124)(3) and 
ms = 2, then a^ = (1524) (3). As described in [H], one can define on Eoo the so-called Ewens 
measure, which is the equivalent on virtual permutations of the Ewens measure on Sat. More 
precisely, let 9 G M+ be a parameter, and let {Mj)j>i be a sequence of independent random 
variables, Mj G {1, ..., j} such that for all j > 2: 



and 

nM, = *1 = ;^ 

(0) 

for all k < j. On Sqo, the Ewens measure /loo of parameter 9 is defined as the image of the 
law of {Mj)j>i by the map $oo- The name of this measure is consistent, since the image of 
yUoo by vToo.AT is precisely the Ewens measure /i}^ of parameter 9 on S^r, under which the 
probability of a permutation a^ is given by the expression: 

^c(o-iv)-l 



i9 + l){9 + 2)...{9 + N-l)' 

where c{a]\f) is the number of cycles of a^- 

Remark: For 9 = 1, the Ewens measure on Sn is the uniform measure, and for ^ = 0, it 
is the uniform measure on permutations with a unique cycle. 

In [11] it is proved that if (crAr)Ar>i is a virtual permutation following the Ewens measure of 
parameter 9, then for all k > 0, the k-th length of cycle (by decreasing order) corresponding 
to the permutation a^, divided by A^, tends a.s. to a random variable Xk when A^ goes to 
infinity. Moreover, the decreasing sequence {xk)k>i follows the Poisson-Dirichlet distribution 
of parameter 9. This property can in fact be easily explained by the construction of the Ewens 
measure on S 00 we give below. 

Let A = {Xj)j>i be a decreasing sequence in R+ and let us denote: 

K{X) := mi{k > 1, Afc = 0} G N* U {cx)}. 

The set E{X) is defined as the disjoint union of circles {Cj)i<j<Kix), such that Cj has perime- 
ter Xj. Now let X = (xfc)fc>i be a sequence of distincts points in E{X). One defines a virtual 
permutation o"oo(A,a;) = {aN{X,x))j^y^ as follows: for A^ > 1, fc G {!,... ,A^}, there exists 
a unique j such that the point Xk lies on the circle Cj. Let us follow the circle Cj, coun- 
terclockwise, starting from x^'- the image of k by o"7v(A,x) is the index of the first point in 
{xi, ...,xn} we encounter after Xk (for example, if Xk is the only point in Cj and {xi, ..., xjy}, 
then fc is a fixed point of <7iy{X, x), because starting from Xk, we do a full turn of the circle Cj, 



before encountering Xk again). The cycle structure of criv(A, x) is the following: two elements 
k and / in {1, ..., iV} are in the same cycle if and only if Xk and xi lie on the same circle, and 
the order of the elements {fci,...,/cp} in a given cycle corresponds to the counterclockwise 
order of the points x^^, ..., x^ , which are on the same circle. Moreover, the cycle structure of 
o"Ar+i(A, x) can be obtained from the structure of a^i^X, x) by a " Chinese restaurant" process: 

• If xn+1 is on a circle which does not contain any of the points xi, ...,XAr, then one 
simply adds the fixed point A^ + 1; 

• If xjsi+i is on a circle which contains some of the points Xi, ...,XAr, and if xat+i lies 
just before Xp if one follows this circle counterclockwise, then one inserts A^ + 1 in 
the cycle containing p, just before p. 

This construction implies that ((TAr(A,x))^>j^ is a consistent sequence of permutations, and 
then croo(A,x) is a virtual permutation. Therefore a virtual permutation can be viewed as 
a Chinese restaurant process with continuous tables (the circles (C'j)i<j<a'{a)), and with an 
infinite number of customers (the points {xk)k>i)'- its component of index A^ is obtained by 
taking into account only the A^ first customers. Note that for the moment. A, the sequence 
of lengths of the circles (or the tables!) plays a minor role in our construction. However, 
it becomes important when one introduces randomness. More precisely, in this framework, 
one obtains the following construction of the Ewens measure on Eqo (and therefore, on E^r, 
by using tToo.n)- 

Proposition 2.1. Let 9 G M+. // the random sequence A follows the Poisson-Dirichlet 
distribution of parameter 6 (for 6 = Q, one sets Ai = 1 and A^ = for k > 1), and if, 
conditionally on \, the points {xk)k>i ore i.i.d., with distribution absolutely continuous with 
respect to the Lebesgue measure (e.g. the uniform distributioiij) on E{\) (and hence, a.s. 
distinct), then for x := {xk)k>i! the virtual permutation a^o (A,a;) follows the Ewens measure 
of parameter 6. 

Proof. The probability law of a random virtual permutation is uniquely determined by its 
image by ^^^ which is the probability law of a sequence (Afj)j>i of random variables {Mj G 
{1, ..., j}). By the monotone class theorem, this law is uniquely determined by the sequence 
of laws of (^j)i<j<Ar, A^ > 1. Now, by applying $Ar for all A^ > 1, one deduces that the 
law of a virtual permutation is uniquely determined by the sequence of laws of its images 
by vTocAT, A^ > 1. This property implies Proposition 12.11 if one shows that for all A^ > 1, 
(Jn (A, x) follows the Ewens measure of parameter 9. To prove this, let us first observe 
that for all permutations ip G Sn, the law of {x^(^k))i<k<N is equal to the law of {xk)i<k<N- 
Hence, the law of a^ (A, x) is invariant by conjugation. This imphes that the probabihty 
of a given permutation, under this law, depends only on its cycle structure, as under the 
Ewens measure. Therefore, it is sufficient to prove that for all partitions (/i, ..., /p) of A^, the 
probability that the supports of the cycles of utv (A, x) are exactly the sets of the form 

{/l + ...+/„ + 1, /l + ... + /„ + 2, ..., /i + /2 + ... + /n + In+l] 



Here, uniform on E{\) means the following: Xk lies in Cj with probability Aj, the perimeter of Cj (note 
that the sum of the perimeters is a.s. equal to one), and conditionally on Xk G Cj, Xk is uniform on the 



circle Cj. 



for < n < p — 1, is the same as under the Ewens measure. Now, conditionally on A, this 
probability can be written as follows: 

Hence one only needs to prove the equality: 



E 






QP-I 



+ i)...{e + N- 



^n(^.-i)! 



(2.1) 



n=l 



since its right-hand side is the probability of the event described above, under the Ewens 
measure of parameter 6. For ^ = 0, both sides are equal to one if p = 1 and /i = N, and 
to zero otherwise: therefore one can assume that ^ > 0. Let {lJ-j)j>i be a random sequence 
following the GEM law of parameter 6, and let {i^j)j>i be the sequence obtained by putting 
(/ij)j>i in decreasing order. One obviously has 



E n^fc= E n 



,^n 



Now, {i^j)j>i and {\j)j>i have the same law, hence it is sufficient to prove that 



E E n^S: 

One needs the following lemma: 



+ l)...{0 + N-l)l\ 



U^in-iy. 



Lemma 2.2. Let {^j)j>i he a GEM process of parameter 9 > Q, and r,s E 
quantity 



(2.2) 



Then the 



E{e,r,s) :=E 






j=i 



satisfies the equality 



Proof: One can write 



E(e r s) = "'(^ + ^-^)'^ 
^ ' '^ (r + s + e-l)\(r + s) 



E{9,r,s)=E[fi[{l-fi, 



+ E 
Now conditionally on /ii, ( j 

E 



l-fiiY+'E 



H 






/ii 



'Ml 



i>i 



. i=2 \ "^ '' \ j=2 

is a GEM process of parameter 6. Therefore: 



fi(r^y(-sA) 



/il 



E(^,r,s). 



which imphes 

E{9,r,s)=E[fi[{l-fi,y]+E{9,r,s)E[{l-fi,r+% 

and, since the density of the law of /ii, with respect to the Lebesgue measure is 6^ (1 — x)^^^ 
on (0,1), 



E{e,r,s) 



EK(1 - fiiY] _ 9 f^ a;*^!! - xY+^-^dx 



1 - E[(l - ^iiY+^] i-ej^{l- xY+'+'-'dx 

er\{s + e- l)!/(r + s + ^)! er\{s + 6- l)!(r + s + 



l-9/{r + s + 



[r + s){r + s + 



which imphes Lemma [2.21 

Now let us go back to the proof of Proposition 12 .![ For all integers g > 0, 1 < ji < ... < jg, 

and ri, ...,rg,r, s > 0: 



E 



E H/^'U^ 



^Jp o 



l-J^/ii 



i=i9+i \p=i 



j=i 



(^i) 



X E 



..?, r 



fij 



j=jq + i- 



Efti Ai 



'-_?,! 



l<i<i. 



/Uj 



n 

s 



fl 



Ip 



3q 
1=1 



r+s 



'«=iq + l 



- Eili /i 



(^ 



i)l<i<jq 



Since conditionally on (yU 



j/i<j<jp) 



/^iq+i' 



-'- Z^i=l l^i J jiyi 



is a GEM process of parameter 9, the last conditional expectation is equal to E{9^ r, s). One 
deduces that 



E 



E n^2u^ 



"Jp o 



1-5^/ii 



i=i<?+i \p=i 



i=l 



E(^,r,s)E 



J9 



n"?.' i-E 



/Uj 



^p=i 



i=l 



By considering all the possible values of (jp)i<p<g, and by adding the equalities, one obtains 

that 



E 



Jq+l 



jl<J2<---<jq+l \P=1 



« = 1 



E{9,r,s)E 



3q 



r+s 



E n^SMi-E". 

jl<J2<---<jq \p=l 



i=l 
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By applying recursively this equality, and using Lemma 12. 2[ one deduces that 



E 



jl<J2<...<jq \P=1 



■Jp 



m=p+l 



n 

p=i 



J^ e {r,)\ \[YL^^^. 



,.fn.A ^^1^ n^ 



We can now compute the left-hand side of (12. 2p : 



E 






Ee 

o-eSp 



E n^ 

<---<V(p)"= 

E ft"' 



(T(n) 



V(l)<---<V(p) " = 1 
P 






E^MH^'^w' 



-11 



1 

11^-1 r« + iv-i)! A.nLiE:.„^. 



n=l / V- ■ -■ -/- ^g£^ xxn=i ^^m=n--M 



Therefore, (12. 2p and then Proposition 12. ![ is proved if one checks the equality: 

1 1 



E 



1 ln=l 2^m=n 'o-(m) 1 ln=l '" 



Now, since for all / > 0, 



one deduces that 

E 



e dx. 



(tSSb 



lln=l l^m=n^<^{'rn) 



E 



" zZn^l ^ri [jZTn^n ^ 



a{m) 



JL~\_ ^ » » » ^JL T}-;^ y) 



JJdx„ 



n=l 



\ / g~ Em = l 'CT(m) En = l ^'" I I j^^; 



By doing the change of variable 



yo{m) = Yl 



Xr. 



n=l 
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one obtains: 

1 r p ^ 

— / p~ 22m = l ^rnVm I I ^„, — 

- / e ll"y" ~ HP / 

which completes the proof of Proposition 12.11 D 

As an illustration, we quickly show how this proposition implies the following almost sure 
convergence result for relative cycle lengths due to Tsilevich: 

Proposition 2.3 (Tsilevich [S]). Let {(yN)N>i be a virtual permutation following the Ewens 
probability measure with parameter 9. One defines the sequence (a^ )k>i of normalized 
lengths of cycles of ajq (i. e. lengths divided by N), ordered by increasing smallest elements, 
and completed by zeros. Then, for all k > 1, a]^ converges almost surely to a random 
variable cr^ , and {a^ )k>i follows a GEM distribution of parameter 6. In particular, the 
law of (a^. )k>i converges weakly to the GEM(6) distribution. Moreover, if y\ denotes the 
l-th largest element of (a^ )a:>i for all integers N >1 and N = oo, then 

(N) ^ (oo) 

Vi ^ Vi 

W— 5>oo 

a.s., and (Vi )i>i follows a Poisson-Dirichlet distribution of parameter 9 . In particular, the 

law of {yl )i>i (i.e. the sequence of decreasing normalized lengths of cycles) tends to the 
PD(9) distribution. 

Proof. Let us construct (crAr)Af>i via Proposition 12.11 Since Proposition 12.31 is trivial for 
= 0, one can suppose 6' > 0. This implies (with the notation of Proposition 12. ip that 
\j > for all j > 1, and a.s., there exists p such that Xp lies on the circle Cj. We define a 
sequence (j„)n>i by the following recursive construction: 

• The index ji is given by: Xi G Cj-^ : 

• For n > 1, ji,...,jn already defined, jn+i is given by: Xp G Cj^^-^, where p is the 
smallest index such that Xp ^ Cj-^ U ... U Cj^^ (this index a.s. exists). 

It is easy to check that for fc > 1: 

(JV) _ \Cj^n{xi,...,XN}\ 1 >^ 

^ ~ N ~ iV ^ "''^^^■^ ■ 

p=l 

Now, by the law of large numbers, it is almost sure that for all integers j > 1: 

1 ^ 

p=l 

Then, a^ tends almost surely to Aj^.. Now, by construction of the sequence {jk)k>i, we see 
that (Ajj.)fc>i is the classical size-biased reordering of (Aj)j>i, and hence a GEM process of 
parameter 6. Now, it is obvious that {yi )i>i is the decreasing reordering of (Aj^)fc>i, i.e. 
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the PD(9) process {Xi)i>i. It remains to prove that y| — )■ A^ almost surely. Indeed, if / is 
fixed, there exists a.s. a (random) index A^o > ^ such that: 

No 

fc=i 
Since a^. — > a^ for all k > 1, there exists a.s. a (random) Ni such that for A^ > Ni. 

(N) 



fc=l 



a^' > 1 - A,/2. 



Since the numbers (a^°° )a;>i are a.s. pairwise distinct, there exists a.s. N2 such that if 
N > N2, the order of (a^ )i<k<No is the same as the order of (a^ )i<A;<Aro- I^i particular 
(recall that Nq > I), the l-th largest element of (a^ )i<fe<Afo has the same index r as the /-th 
largest element of (a)^ )i<k<No- Now, since 0;^°° < Xi/2 for all k > Nq, ar is also the l-th. 
largest element of (ctj^ )k>i, i-e. A^. For A^ > sup(A''i, A^2)5 ctr is the l-th largest element of 
{al )i<k<No and a\. < Xi/2 for A^ > A'o; hence the l-th largest element y} of {a}, )k>i 



is included in the interval [ar , Or V A//2]. Since ar — > ar = Xi, Proposition 12.31 is 



.(^\«(^)vA,/2]. Since a(^) -^ ai^^ 

N^oo 

proved. D 

3. The "non-unitary case" 

3.1. The normalized and non-normalized empirical eigenvalues distributions. Let 

M be a matrix in ^(A^) for some A^ > 1. We associate with the point process of the 
eigenvalues of M the finite measure fi{M) on C defined by 

l^iM) := Yl ^^MiX) 6x, 

where E(M) is the set of eigenvalues of M , mM{X) is the multiplicity of A as an eigenvalue of 
M , and 5\ is Dirac measure at A. By the general description of eigenvalues given in Section 
II. 2^ one has, for all a G S^r, Zi, ..., ^at G C*: 



li{M{a,zi,...,ZN)) = Y 5Z ^' 



771=1 ijjlr. 



where /i, I2, ■■.,ln are the lengths of the cycles Ci, C2..., C„ of a, and for 1 < m < n: 

jeCm 

Let us now suppose that the distribution of a sequence of random matrices {Mn)n>i, Mn G 
^(A^), is of the form P(oo,^,£). One has Mj<i = M{a]\[,zi, ...,zn) where (o"7v)Ar>i follows 
the Ewens(6') distribution, and is independent of the sequence {zj)j>i of i.i.d. variables, 
which have law C. Since for all A^ > 1, the cycle structure of a^ can be deduced from 
the cycle structure of o"iv+i by removing A^ + 1, there exists a partition 11 of N* such that 
for all A^, the supports of the cycles of a^ are obtained by intersecting the sets of 11 with 
{1,...A^}. Moreover, under the Ewens(^) measure, 11 contains a.s. an infinite number of 
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sets (see Section [2]): let us order them by increasing smallest elements, and denote them by 
{Cm)m>i- One then has 






(3.1) 



m=l 

where lN,m is the cardinality of C^^rn, the intersection of Cm and {1, ..., A^}, and 

The natural question one can now ask is the behaviour of the measure ^{Mj^) for large N . 
Since each cycle of aN gives a number of eigenvalues equal to its length, one can expect that 
^{Mn) is dominated by the large cycles of M^- Moreover, the / eigenvalues corresponding 
to a cycle of length / form a regular polygon of order /, and the distance of their vertices to 
the origin is equal the the l-th root of the product of / i.i.d. random variables of law C. If 
/ is large and if one can apply a multiplicative version of the law of large numbers, one can 
expect that this distance does not vary too much. Then, it is natural to guess that under 
some well-chosen conditions on C, the measure //(M^v), which has total mass A^, is close to 
A^ times the uniform measure on a circle centered at the origin. Indeed, we can prove the 
following statement: 

Proposition 3.1. Let {Miy)N>i be a sequence of matrices following the law F{oo,6,C) for 
some 9 > and some probability C on C* . We suppose that if Z is a random variable which 
follows the distribution C, then log(|Z|) is integrable. Under these assumptions, almost 
surely, the probability measure n{M]sr)/N converges weakly to the uniform distribution on 
the circle of center zero and radius exp (E[log(|Z|)]). 

Proof. Let / be a continuous and bounded function from C to M, and let i? > 0. Then, there 
exists a constant A > 0, and a function a from (0, R) to M+, tending to zero at zero, such 
that for all e G (0, R), for all integers / > 1, and for all z G C such that \z\^^^ E {R — e, R + e): 

r.27r 



E 



/H 



2^ 



f{Re'^)dX 



< - + la{e) 



(3.2) 



Indeed, let us define for all 5 > 0: 

r^{5) := sup{|/(y) - f{y% \y - y'\ < 5, |y|, \y'\ < 2R}, 

which tends to zero with 6 since / is uniformly continuous on any compact set. With this 
definition, we obtain: 



E /(^) - E /H 



< ivi^) 



where z' := zR^ /\z\ has modulus RK Now, there exists A G [0,27r//) such that: 

5^/(a;) = ^/(i?u;e^^)=:$(A) 



a;'=l 



One has, for all A, A' G [0, 2tx/1): 



|<I>(A) - $(A')| < lr]{R\\ - X'\) < lr]{2'KR/l). 
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Moreover 



2tt/1 i>2n 

$(A)rfA= / f{Re'^)d\ 
Jo 



and then, for all A G [0, 271/1), 

HX) - ^ I ^ f{Re'')dX 



< lri{2nR/l), 



which implies: 






< / [r^(e) + r/(27ri?//)] 



If / < 27ri?/e, one can majorize this quantity by 47ri??7(27ri?)/e, and if / > 27ri?/e, one can 
majorize it by 2lri{e). Hence we obtain (13.21) . 

Since for S > depending only on /, the left-hand side of (13. 2p can be trivially majorized 
by Bl for any 2; G C, we deduce, for /i equal to the uniform measure on the circle of radius 
R: 



l^jfdfxiM^)-Jfdfx 



WjV.m>0 



m=l 



Bl 



N.m 



1, 



|i/ijv 



AT |^iv,mr"^'™^{fl-€,K+e) ' Ne N 



^ , ''N,m , N 



For now, let us take R := exp (E[log(|Z|)]), where Z is a random variable following the law 
C By the strong law of large numbers applied to the sequence (log |^j|)jgC™, it is not difficult 
to check that a.s., for all ?7i > 1: 



Bl 



N,m 



1, 



tends to zero when A^ goes to infinity. Moreover, independently of A^, this quantity is 
dominated by Bsm, where Sm is the supremum of lN,m/N for A^ > 1. For the moment, let 
us assume that a.s.: 



Sm < 00 



m=l 

In this case, one can apply dominated convergence and obtain: 



(3.3) 



m=l 



Bl 



N.m 



N,m>0,\Zr,,^\'-^''^'"^t(R~e,R+e) jy 



0. 



Af^oo 



Moreover, 



m=l 



A 

' Ne N~^oc 



E ^'iv,m>0^r^ ^ 



a.s., since the number of cycles of a^ increases slowly with respect to N (the order of 

magnitude is log(A^)), and 

00 , 

E^ ''N,m 

^lN,m>0' 

m=l 

Hence we deduce that 



N 



a[e] = aie). 



lim sup 

Af-s>oo 



ljfdf,iM,,)-jfdfi 



<«(e), 
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and by taking e — )■ 0, we are done. It only remains to prove (13.31) . This relation can be shown 
by looking carefully at the construction of the Ewens(6') measure on virtual permutations 
given in Section[2l Indeed, for A^ > 1, conditionally on {lK,m)m>i,i<K<N , with mo := inf {m > 
1, ^7v,m = 0}, one has Iw+i^m = ^7v,m + lm=mi whcre mi is a random index equal to m' with 
probability lN,m/{N + 6) for m' < mo and to rriQ with probability 6/{N + 6). This implies 
quite easily that {lN,m/{N + 6))n>i is a nonnegative submartingale. By Doob's inequality 
one deduces that the expectation of s^, is dominated by a constant (depending only on 0) 
times the expectation of y^, where Hm is the limit of lN,m/N, which exists almost surely. 
Now, {ym)m>o is a GEM process of parameter 9, hence the expectation of y'^ decreases 
exponentially with m. Consequently, 



E 

m=l 



E[s„] < oo, (3.4) 



which implies (13. 3p almost surely. D 

If we do not want to deal with virtual permutations, we can replace the a.s. convergence by 
a weak convergence in probability, as follows: 

Corollary 3.2. Let {Mn)n>i be a sequence of matrices such that Mjv G ^(A^) follows the 
distribution F{N,6,C) for some 6 > and some probability C on C* . We suppose that if 
Z is a random variable which follows the distribution C, then log(|Z|) is integrable. Under 
these assumptions, the probability measure fi{M]\f)/N converges weakly in probability to the 
uniform distribution on the circle of center zero and radius R := exp (E[log(|Z'|)]), i.e. for 
all continuous, bounded functions f from C to M.- 

f.27r 



TV Jc ^^°° ^^ Jo 



in probability. 



This convergence result means that most of the eigenvalues of a matrix Mn following 
F{N,9,C) are concentrated around the circle of radius exp (E[log(|Z|)]). Now, even for N 
large, there remain some eigenvalues which are far from this circle, since the law of large 
numbers involved in the proof of Proposition 13.11 does not apply for the small cycles of the 
permutation a^ associated with Mj^. In order to study the influence of the small cycles, let 
us suppose that (Mjv)Ar>i follows the measure P(cxd, 9, C) (which is possible since the image 
of P(oo, 6, C) by the A^-th coordinate is P(A^, 6, £)). Then, one can write the measure /i(Mjv) 
in the following way: 

oo 

^l{MM) = Y, Y. Y. -^^ 

fc=l m>l,i]V,,„ = fc(^fc=Z]v,m 

with the same notation as in equation (13.11) . This equality implies the following equality in 
distribution: 

oo O-N.k 

f,{M^) = Y,Y. E '^- (3.5) 

k=l p=l <^'==Tfc,p 

where for all A; > 1, a7v,fc is the number of A;-cycles of the permutation cttv (which follows the 
Ewens(^) measure on S^r), where for k,p > 1, the law of T^p is the multiplicative convolution 
of k copies of the distribution C, and where {aN,k)k>i and the variables T^p, k,p > 1 are 

16 ' ~ 



independent. Now, the finite dimensional marginals of {aN,k)k>i converge, in distribution, 
to the corresponding marginales of {ak)k>i, where the variables a^ are independent Poisson 
random variables, with E[afc] = 6/k (see for instance [2]). One can then expect that in a 
sense which needs to be made precise, the law of fi{M]\f) converges to the distribution of 



/ic 



oo afc 

fe=i p=i 1 



E ^^ 



(3.6) 



k>p 



where all the variables a^ and Tf^^p in sight are independent. Of course, one needs to be 
carefull, because the measure fi^o has an infinite total mass, which, under the assumptions of 
Proposition l3.lt is expected to concentrate around the circle of radius exp (E[log(|Z|)]). One 
also remarks that the convergence expected here is very different from the convergence proved 
in Proposition 13. Ij in particular, it involves the measure ^{Mn) and not the probability 
IJ,{M]S[)/N . In order to state a rigorously our result, let us give the following definition: 

Definition 3.3. Let X be a real, integrable, random variable. For g > 0, we say that X is in 
W if and only if for a sequence {Xk)k>i of i.i.d. random variables with the same distribution 
as X and for all e > 0, there exists C > such that for all ra > 1: 



P 



^E^^)-E[x] 



fc=i 



> e 



C 
< — . 



Remark 3.4. If g > 1 is an integer, by expanding 



E 



\ k=l 



29 



E[X]) 



and by using Markov's inequality, one easily proves that a random variable in L^"^ is also in 
W. 

We can now state a result about the convergence of ^{M^)- 

Proposition 3.5. Let {Mi\f)N>i be a sequence of matrices such that Mn G Q{N) follows the 
distribution ¥{N, 6, C) for some 6 > and some probability C onC* . We suppose that if Z is 
a random variable following the distribution L, then log(|Z|) is in W^ for some g > 0. Under 
these assumptions, for all bounded, continuous functions f from C to M, such that / = 
on a neighborhood of the circle \z\ = R, where R := exp (E[log(|Z|)]), / is a.s. integrable 
with respect to fi^o, where the random measure /ioo is given by f l3.6p . Moreover, the following 
convergence in distribution holds for such test functions f: 



fdfi{M, 



N 



N-^oo 



fdlXc 



Proof. Let (^r)r>i be a sequence of independent Bernouilli random variables, such that the 
parameter of ^r is equal to 9/ {9 + r — 1) (in particular, ^i = 1 almost surely). We suppose 
that {C,r)r>i is independent of {Tk^p)k,p>i, and for all N,k > 1, we define bN,k as the number of 
pairs of consecutive ones in the sequence (^1,^2, ■■■,^n, 1) which are separated by a distance 
of k, and bk as the analog for the infinite sequence {^r)r>i- In other words: 



^N,k 



1 



N+l~k>l,(N+l_k = l, 5iV+2-fc = 



+ 



E 



1 



§j=5i+fe=l'?j+i = 



-ij+k- 



^<3<N-k 
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and 

By the classical properties of the Feller coupling (see for instance [2j), {bk)k>i has the same 
distribution as {ak)k>i and for all N > 1, {bN,k)k>i has the same distribution as (aAr,fc)fc>i- 
Therefore, in Proposition 13. 5[ one can replace fi^M^) by un and /Xoo by z/qo, where: 

oo t'jv,fe 
k=l p=l a;fe=Tfe,j, 

and 

oo bk 

^-•=Z1Z1 Yl '^- 

fc=l p=l a;fc=Tfc,p 

Let / be a continuous and bounded function from C to M, equal to zero in a neighborhood 
of the circle \z\ = R. There exists A > 0, < Ri < R < R2, depending only on /, such that 

\f{z)\ < A1|^|^(R,,R2), and then, 

00 6, 



„ 00 Ofc 

/ \f\du^ < AJ2J2''hn.r,\'/''m■uR2)■ 



fc=i p=i 

The function / is a.s. integrable with repect to z/qo (and then, /ioo)- Indeed, 

/ \f\di^oo < 00 
Jc 

if and only if 

00 bfc 

fc=l p=l 
and the expectation of this quantity is 



"E 



A; 



fc=i 

which is finite since log(|Z|) is in W for some g > 0. Let us now introduce the random 
measure: 

""^ = ^2^ Yl '^-' 

fc=l p=l a;'==Tfc,p 

where 

CN,k = 2^ l5^.=5^.+^=i,g^.+i=...=5^.+^_i=o- (3.7) 

One has feAr^fc = CN,k, except for k equal to the smallest integer fco such that ^N+i-ko = 1, in 
which case b^^ko = CN,ko + 1- Therefore, 
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where 

Since / is a.s. integrable with respect to Uoo, and for all k > 1, bk is the increasing limit of 
cjvfe, one obtains, a.s.: 



fdu', 



N 



N^oo 



f duo 



{3.t 



Moreover, there exist < -Ri < -R < R2, depending only on /, such that for all k' > 1: 



P 



fdu'1,^0 



k>l 

< P[fco < k'] + sup P f|T,»,i|i/^" ^ (i?i,i?2) 



k">k' 



Now, P[A;o < k'] tends to zero when A^ goes to infinity (the order of magnitude is 1/A^), and 
by taking N ^ 00, and then fc' — )■ 00, one deduces: 



P 



fdu';,^o 



0. 



Therefore, by taking (13. 8 p and (13. 9p together: 



fdi^'r 



N 



7V->oo 



Af-s>oo 



fduo 



(3.9) 



in probability, and a fortiori, in distribution. 



D 



Remark 3.6. In the proof of Proposition 13. 5[ we have used the Feller coupling in order to 
replace the convergence in law by a convergence in probability. However, this coupling does 
not correspond to the coupling used by considering the measure F[oo,9,£). Moreover, if 
(MAr)7v>i follows this measure, the number of cycles of a given length in a^ does not converge 
when N goes infinity, and the support of these cycles is changing infinitely often. Hence, 
one cannot expect an almost sure convergence (or even a convergence in probability) in 
Proposition 13. 5[ 

Propositions 13.11 and 13.51 apply for a large family of distributions C, however, some integra- 
bility conditions need to be satisfied. One can ask what happens if these conditions do not 
hold. We are not able to prove a result for all the possible distributions C, but it is possible 
to study some important particular cases, if the probability distributions involved in our 
problem can be explicitly computed. Here, the most important distributions we need to deal 
with are the probability measures Ck, k > 1, defined in the following way: Ck is the unique 
measure on C* , invariant by multiplication by a k-th root of unity, and such that its image 
by the k-th power is the multiplicative convolution of k copies of C Intuitively, Ck is the law 
of a random A;-th root of the product of k independent random variables with law C, chosen 
uniformly among the k possible roots. The measures Ck are not easy to compute in general. 
One case where the computation is simplified is the case where C has the radial symmetry. 
Indeed let (l^)p>i be a sequence of i.i.d. random variables with the same law as log(|Z|), 
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where Z follows the distribution C The distribution L^ is the law of e*®'*"^^p=i^'' where 6 
is independent of (^)p>i and uniform on [0,27r). In particular, if L is the law of e*®+^'^'*, 
where p is a strictly positive parameter, B a uniform variable on [0, 27r) and Sa [pi € (0, 2]) 
an independent standard symmetric stable random variable of index a, then L^ is the law 
of e^^^P^ ° °'^". Using this explicit description of £fc, we can make a detailled study of 
the "stable case". For a > 1 (and in particular for a log-normal modulus, corresponding to 
a = 2), Propositions 13 . 1 1 and 13.51 directly apply. Therefore, let us suppose a < 1. For a = 1, 
one has the following: 

Proposition 3.7. Let p > 0, and let C be the law o/e*®^'''^^ where Q is a uniform random 
variable on [0, 27r) and Si an independent standard symmetric Cauchy random variable. For 
6 > 0, let {Mn)n>i be a sequence of random matrices such that Mn follows the distribution 
¥{N,6,C). Then the distribution of the random probability measure p{Mn)/N converges to 
the law of the random measure 



Eu 



(m) ) 

m>l 

where {xm)m>i is a Poisson-Dirichlet process with parameter 9, (S*}™ )m>i is an independent 
sequence of i.i.d. standard symmetric Cauchy variables, and for R> 0, fi^ is the uniform 
measure on the circle of center zero and radius R. This convergence has to be understood as 
follows: for all continuous and bounded functions f from C to M, 

fdfi{MN) — y fdfi^ 

Jv Jc 

in distribution. 



N Jq Af-s>oo 



Proof One can suppose that (M7v)Af>i follows the distribution P(oo,6',£). With the same 
notation as in the proof of Proposition 13. ![ one sees that p{Mm)/N has the same distribution 
as the random measure: 



iNm. I 1^ 



n^ -f:. — 2^ '^-^ 



LdLirn 



\ ''N.m 



where {Lm)m>i is a sequence of i.i.d. random variables with law C (recall that Ck is equal 
to C for all A; > 1). If / is a continuous, bounded function from C to R, one has: 



'^'^ m=l \ 



l-N.m ^j, 

Now, there exists a GEM process {ym)m>i of parameter 6 and a sequence of random variables 
(sm)m>i, such that almost surely, 

oo 

Sm < OO 



m=l 



'N,m 

AT- ~ ^ ymi 
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for all r/i > 1 and 

l'N,m ^ 

for all A^, -m > 1 . By the convergence of Riemann sums corresponding to the integral of 
continuous functions on a compact set, one deduces that almost surely, for all m > 1: 



N \ lN,m ^„, 

'■N,m 



.f-^ / iv^oo 27r Jo 



where the left-hand side is smaller than or equal to Sm||/||oo; independently of A^. By 
dominated convergence: 

f.27r 



If N^Qo Zn ^ — ' In 



e'^)dX, 



almost surely. This implies Proposition 13.71 D 

Remark 3.8. Almost surely, the random measure /ioo is strictly positive for any nonempty 
set of C Therefore, for all continuous functions / from C to M+, non-identically zero, and 
for all A G M+: 



P 
Hence, one cannot expect an ana^ 



f fdfi{MN)<A 

Jc 



— > 0. 



og of Proposition 13.51 in the case studied here. 



For the case a < 1, one intuitively expects that most of the eigenvalues become very large 
or very small. The precise statement is the following: 

Proposition 3.9. Let p > 0, and let C be the law of e'^®^^^°' , where Q is a uniform random 
variable on [0, 2tt) and Sa an independent standard symmetric stable random variable, with 
index a < 1. For 9 > 0, let {Mn)n>i be a sequence of random matrices such that Mj^ 
follows the distribution P(A^, 6, C) . Then, the distibution of the random probability ^{Mjq)/N 
converges to the law of the random measure GqSq, where Gq is a beta random variable with 
parameters {9/2,9/2), in the following sense: for all continuous functions f from C to M, 
with compact support, 

1 [ fdf,{M^) -^ GefiO) 

N Jc Af->oo 

in distribution. 

Remark 3.10. The total mass of the limit measure is a.s. strictly smaller than one. Intuitively, 
this is due to the fact that a large part of the total mass of ^{M^) is going far away from 
zero, when A^ is large. This mass is missing in the limit measure in Proposition 13. 9[ because 
we consider functions / with compact support. 

Proof. We suppose that {M]s[)n>i follows the distribution P(oo, 9, C). Let / be a continuous 
function from C to M, with compact support, and let us choose R > 1 such that f{z) = {) 
for all z such that \z\ > R. Let us define, for all 2; G C, 



and for r G (0, 1): 



(?(2;) := l|,|<i/(0). 

/3(r):=sup{|/(^)-/(0)|,|^|<r}, 
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which tends to zero with r. One checks that 



N 
which imphes: 



\f-g\dfi{M,,)</3{r) + 2\ 



/i(M^)({^e C,r < \z\ < R}) 

' N ' 



P 



jj J\f-9\ d^iiM^) > m + (2II/II00 + l)r 



< P 



r 



fi{MN){{z EC,r < \z\ <R}) 
iV 
fi{MN){{ze C,r < \z\ < R}) 

N 



> r 



Now, 



E 



7i(Mjv)({^eC,r < \z\ <R}) 

iV 



m=l 



N,m 



N lN,m>0,\ZM.m\'-''^''^-"'^&{r,R) 



From the independence of 1^^^ and {zj)j>i, and from the basic properties of stable random 
variables, one has 

P [lN,m > 0, IZiv.ml^'/'^''"^ e (r,i?)|/jv,m] = ^(/^.m), 

where the function \E' (which can depend on a, p, r and R, but not on A^) is bounded by one 
and tends to zero at infinity. Hence: 



E 



)i{MN){{z e C,r < |;z| < R}) 



N 



< E 






lN,m^{lN,m) 

N 



Now, for all m > 1, a.s.: 



and 



By dominated convergence, 



E 



lN,m^{l 



N,m) 



N iV-5-oo 



A^ 



< Sr 



fx{MN){{zeC,r< \z\ <R}) 

N 



N^oo 



which implies: 



P 



A^ 



|/-(7|d/i(Miv)>/3(r) + (2||/|U + l)r 



N-^00 



0. 



By letting r go to zero, one deduces that 

1 



in probability. Therefore, it is sufficient to prove the conclusion of Proposition 13.9^ with / 
replaced by g. Moreover, one can suppose /(O) = g{0) = 1. In this case, one has: 



^^ m>l 



'■N,m ^ 

AT ^]V,m>0,|^Ar,ml<l' 
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Now, by symmetry of the stable variables considered here and by the independence of 

{^N,m)m>l and {Zj)j>l, 



in distribution, where (em,)m>i is a sequence of i.i.d Bernoulli variables of parameter 1/2, 
independent of {lN,m)N,ra>\- Now, by dominated convergence (recall that the sum of the 
variables Sm is a.s. finite), 

N N-^oo ^-^ 

»n>l m,>l 

which implies Proposition 13.91 if we check that 

is equal to Gq in distribution (recall that {ym)m>i is a GEM(6') process, independent of the 
sequence (em)m>i)- This fact can be proved as follows: by self-similarity of the GEM(6') 
process, 

X = VX + e{l-V), 

in distribution, where V, X, e are independent, l^ is a beta variable of parameters 6 and 
1, and e is a Bernoulli variable of parameter 1/2. With this identity, one can compute the 
moments of X by induction, and finally, one can identify its law. D 

Again for the stable case with parameter a < 1, one has an analog of Proposition 13.51 for 
the eiganvalues which are in a compact set not containing zero. The precise statement is: 

Proposition 3.11. Let p > 0, and let C be the law of e^^^^^" , where Q is a uniform random 
variable on [0, 2tt) and Sa an independent standard symmetric stable random variable, with 
index a < 1. For 6 > 0, let (Mjv)jv>i be a sequence of random matrices such that Mjv follows 
the distribution ¥{N, 6, C). Then, for all continuous functions f from C to M, with compact 
support, and such that f = on a neighborhood of zero, f is a.s. integrable with respect to 
Hoo, where the random measure /ioo is given by (13. 6p . Moreover, the following convergence 
in distribution holds with such test functions: 



j fdp{Mr,) —^ f fdfi, 

Jc ^^°^ Jc 



Proof. The proof is similar to the proof of Proposition 13.51 One only has to change the 
estimate: 

nTk,i\'^''^{Ri,R2)]<Ck-'^ 
(for some C, g > 0), by the estimate: 

nTk,i\'^''e{r,R)]<Ck-'i, 

available for all r, R such that i? > r > 0, by the classical properties of symmetric stable 
random variables. D 
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3.2. The average eigenvalues distributions. Another interesting problem is the study 
of the expectation jin of the random measure //(Mtv), where M^ follows the distribution 
P(A^, 6, £). The measure fiN of a Borel set in C represents the expected number of eigenvalues 
of Mat (with multiplicity) contained in this set, and it can be explicitly computed: 

Proposition 3.12. Let N > 1, 6 > 0, C a probability measure on C* . If M^ follows the 
distribution F{N,9,C), then the expectation fiN of fi^M^) can be represented as follows: 



/iAT 



N 



N{N-l)...{N-k + l) 



C,. 



Proof. Let / be a Borel function from 
one has: 



^ (A^ - 1 + e){N - 2 + e)...{N -k + 

to M+. By taking the same notation as in fl3.5p . 



Jc 



fdfi 



N 



E 



EE E /M 

k=l p=l ujk=Tk^p 



Now, for all k,p>l: 



E 



E /H 



l-k.p 



k I f dCk 



Moreover, (Tk,p)k,p>i is independent of (av,fc)A:>i- One deduces that 

oo 0,N,k 



E 



EE E /HiK^ 

fc=l p=l a;'==Tfe,j, 
oo „ 

^-1 -'C 



k)k>l 



Therefore 



Hn 



^A;E[aAr_fc]£fc. 



fc=i 



By doing explicit computations of E[aAr^A;] (see e.g. p]), one deduces Proposition 13. 12[ D 

Remark 3.13. When the distribution of the Zj^s is the Dirac mass at 1, then Ck is the uniform 
measure on the k-th roots of unity. 

One has a similar result for the limit random measure fioo of ^{Mn)\ 

Proposition 3.14. Let 6 > 0, C a probability measure on C* , {Mj\f)]^>i a sequence of 
random matrices such that M^ follows the distribution F{N,6,C). Then the expectation of 
the random measure fioo, given by (13.61) (it can be considered, in some sense, as the limit of 
Ii{Mn) for N — )• oo^, is the measure: 



Aic 



"E^ 

fc=i 



k- 



Proof. The proof of Proposition 13.141 is exactly similar to the proof of Proposition I3.12[ D 
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Now, we have a sequence of finite measures JIn, defined as the expectation of ^{Mn) 
and explicitly described, and an infinite measure /too, defined as the expectation of /loo- 
Moreover, we know that, for a large class of probability laws C, the random probability 
measure h{Mn)/N converges weakly in probability to the uniform measure on a circle, and 
in a sense which can be made precise, //(M^v) tends to /ioo- Hence, we can expect analog 
convergences for the sequence of measures {fiN)N>i- One indeed has the following result: 

Proposition 3.15. Let 6 > 0, C a probability measure on C* , (MAr)7v>i a sequence of 
random matrices such that Mn follows the distribution ¥{N,6,C). We suppose that if Z 
follows the distribution C, then log(|Z|) is integrable, and we define R := exp (E[log(|Z|)]). 
Then the probability measure j1n/N (which represents the probability distribution of a random 
eigenvalue of M^, chosen uniformly among the N possible eigenvalues), converges weakly 
to the uniform distribution on the circle {z = R}. Moreover, if log{\Z\) is in U^ , then for 
< Ri < R < R2, the restriction of fiN to the set {\z\ ^ (i?i,i?2)} converges weakly to the 
corresponding restriction of fiao, which is a finite measure. 

Proof. Let / be a continuous and bounded function from C to M. By Proposition 13. II (more 
precisely, by Corollary 13 ■2p : 

i- [ fdfi{MM) -^ ^ [^ f{Re^')dX 

in probability. Since / is uniformly bounded, one deduces the first part of Proposition 13.151 
by taking the expectation. Let us now prove the second part. We now suppose that / is a 
bounded and continuous function from C to M, equal to zero in a neighborhood of the circle 
{|2;| = R}. By Proposition 13. 12[ it is sufficient to prove that 



where 

tN,k 



00 „ 

y2\tN,kAN>k-l\ / \f\dCk -^ 0, 

^-^ /^ Af— s>oo 

fe=l ''^^ 

N{N-l)...{N-k + l) 



{N -i + e){N-2 + e)...{N -k + ey 

Each term of the sum converges to zero when N goes to infinity. Hence by dominated 
convergence we are done if we prove that: 



00 „ 



^J/|rf£fc<oo, (3.10) 

fe=i 

where 

Vk,e = 1 + sup{tN,k,e, N > k}. 
Now, for ^ > 1, Vk^0 = 2, and for 6 < 1: 

k\{9-l)\ 

''''' - ^ + ik-i + e)V 

which implies that v^^e is dominated by k^^~^^+ for fixed 9. On the other hand, the fact that 
log(|Z|) is in U^ implies that the integral of |/| with respect to Ck decreases with k at least 
as fast as 1/fc^. This implies (I3.10p . D 
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As above, it is interesting to see what happens if £ is the distribution of e'®+'''^°, where p > 0, 
is a uniform random variable on [0, 27r) and Sa is an independent standard symmetric 
stable random variable, with index a G (0,2]. For the log-normal case a = 2, Proposition 
13.151 applies directly. For the case a G (1,2), one can apply the first part of Proposition 
I3.15[ which gives the convergence of the average empirical measure fl]\f/N, but the second 
part cannot apply. Indeed, by using the classical tail estimates of stable random variables 
one checks that for all nonempty open sets A G C: 

oo 
k=l 

For a = 1, one has Hk = L for all /c > 1, and then: 

/iAT = N L 

for all A^ > 1. The most interesting case is a < 1. 

Proposition 3.16. Let p > 0, and let C be the law of e^^^^^" , where Q is a uniform random 
variable on [0, 2it) and S^ an independent standard symmetric stable random variable, with 
index a < 1. For 6 > 0, let (Mjv)Ar>i be a sequence of random matrices such that M^- 
follows the distribution F{N,6,C). Then the probability measure Pn/N (which represents 
the probability distribution of a random eigenvalue of Mj^r, chosen uniformly among the N 
possible eigenvalues), converges vaguely to half of the Dirac measure at zero. Moreover, for 
all r,R such that < r < R, the restriction o/ /loo to the set {\z\ G {r,R)} is infinite if 
a > 1/2, finite if a < 1/2, and for a < 1/(2 V (3 — 9)) (in particular for a < l/S), it is the 
weak limit of the corresponding restriction of fiN when N goes to infinity. 

Proof. Let / be a continuous function from C to M, with compact support. By Proposition 
13. 9[ one has the convergence in distribution: 

/d/x(M;v) -^ GefiO). 



Since / is uniformly bounded, one can take the expectation and obtain: 

1 f fdpN -^ ^/(O), 

which gives the first part of Proposition I3.16[ Now, one checks that for all r, R such that 
< r < R, Ck{{z G C, \z\ G (r, R)}) decreases as fc^~^/° when k goes to infinity, which gives 
the condition for the finiteness of /loo restricted to {z G C, \z\ G {r,R)}. The proof of the 
convergence of the restriction of JIn to this set toward the restriction of /ioo is similar to the 
proof of the second part of Proposition 13. 151 One only needs to check that: 

oo 

'^VkM^kiiz G C, \z\ G (r,i?)}) < cx), 
fc=i 

which is true, since Vkfi is dominated by k^^~^^+ and Ck{{z G C, \z\ G {r,R)}) is dominated 
by A;^-i/°, for 1/a > 2V(3-^). D 
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3.3. The g-correlation. The measures JIn and jloo give the average repartition of the eigen- 
values (and its limit for large dimension) of a matrix following the distribution F{N,6,C). 
It is interesting to generalize this study to the g-correlation of eigenvalues for all strictly 
positive integers q, i.e. the distribution of the possible sequences of q eigenvalues. More 
precisely, if Mjv follows the distribution P(A^, 6, C), let us consider, for all q > 1, the random 
measure on C: 



/.M(M^)= Yl ^K.--.,)' 

where (i^j)i<j<Ar is a sequence containing all the eigenvalues of M^, with multiplicity. The 
"g-correlation measure" is defined as the average of //['^^(MAr): it is a finite measure jl^ on 
C, with total mass iV!/(A^ — g)!. One checks that fi^''^{M]\[) has (by taking the same notation 
as in (13. 5p ) the same distribution as the random measure 

^N ■= 2^ Z^ ^ ((-^A;,p)fc>l,l<p<ajv,fe) , 

(«'fe,p)fc>l,l<P<ajv feeH^(g,(aiV,fe)fe>l) (E^:,p&S^:,^,. JTk,p)) 

' V '^ ^ fc>l,l<p<ajY ^ 

where W{q, (aAr,fc)fc>i) is the set of families {'Wk,p)k>i,i<p<aMk ^^ nonnegative integers, with 
total sum q, Sk,wk (^fe,p) is the family of subsets of k-th roots of Tk^p, with cardinality Wk^p, 
and A ((^fc,p)fc>i,i<p<ajv,J is defined by: 



A ((-Efc,p)fc>i,i<p<aiv,fe) - / , ^' 



{^<T{j))l<j<q 

for a sequence {ojj)i<j<q containing all the elements of the sets E^p, each element appear- 
ing a number of times equal to the number of sets E^p containing it. Let us now take 
the conditional expectation with respect to {aN,k)k>i- One obtains the following random 
measure: 

E[un I {aN,k)k>l\ = ^ ^C'"^] ((^fc,p)fc>l,l<P<ajv,fe) , 

(«'fe,p)fe>l,l<P<ajv,fe^^(9'('*iV,fe)fc>l) 

for 

C^"^ ((Wfc,p)fc>l,l<p<ajv,J = ff ^ Yl ^■^^'^^ {i^k,p)k>l,l<p<ar,,k) , 

where 

a.&^ ((^fe,p)fc>l,l<p<ajv,fc) 

is the image, by the permutation a of the coordinates, of 

>CM {{Wk,p)k>l,l<p<a^^,) = (g) 4^""^^, 

A;>l,l<p<ajv,fc 

for 
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where {zi,Z2, ..., Zw^ ).Ck is the image of Ck by the apphcation z i-> {zzi, zz2, ..., zzw,^ )■ By 
defining Afc,r as the number of indices p such that Wk^p = r, one deduces: 



E[um I {aN,k)k>i] 



E 



nfc>i (^n;i 



(,'^fc,rjfe>l,r>0fc-'^(^<?il'*iV,fc;fe>lj 



E- 






.A;>l,r>0 

where L (g, (aAr,fc)fc>i) is the set of famihes {\k,r)k>i,r>Q of nonnegative integers such that 



E --^ 



fc,r — 9 



A;>l,r>0 



and for all A; > 1: 



7 ^ Afc^r — OAT^fc- 
r>0 



Note that for all /c > 1, the measure £|, is trivial (it is the unique probability measure on a 
space with one element), hence, one can remove it in a tensor product. One deduces: 



E[l/Ar I {qn, 



k k>l\ 



E n 



Ovfc! 



(Afc,r)fe,r>ieL(g) ^> 



i(''",i-E..>iV)!n,->,iA^,.K''!)^"i ■■■ 



.^. E- 



o-es„ 






^fe,r 



.fc>l,r>l 



where I/(q') is the set of families {\k,r)k,r>i of nonnegative integers such that 

y^ r\k,r = q 



fc>l,r>l 



and where the inverse of the factorial of a strictly negative integer is considered to be equal 
to zero. Then: 



fJ-N 



E 



u 



[Er>l^fc,r]j^>, 



1 



N,e 



(^k,r)k,r>l&L{q) 



nt.>i |A.,.!(r!)^'-l 



E- 



o-es„ 



(g) (4' 



.k>l,r>l 



where, for all sequences of nonnegative integers {Xk)k>i for which the set of k such that 
Afc > is finite: 



u 






E 



n 



0'N,k- 



.k> 



\ {aN,k — Afc)! 



(recall that this quantity depends on 6, as the law of {aN,k)k>i)- By elementary combinatorial 
arguments, one can prove that 



u 



lx.].>_._UfjF'\N + l-j) 



N,( 



Uk>ik'^ 



■V, 



where V is the probability that in a random permutation following the Ewens(^) distribution, 
the integers from 1 + Yli<k-^j ^^^ Yli<k-^j ^^^ ^^ different k-cjcles, for all k > 1. By the 
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Feller coupling, one deduces that 

l^k]k>l 



u 



N,e 



1 



Efc>ifeA,<JV 



n 



N + 1 



Qj2k>l^k 



N + e-j Uk>ik'^' 



Then by denoting L{N, q) the set of families {Xk,r)k,r>i of nonnegative integers such that 

y^ r\k,r = q 



fc>l,r>l 



and 



one obtains the following result: 



fc>l,r>l 



kXk^r < N, 



Proposition 3.17. Let 6 > and let C be a probability distribution on C*. If one takes the 
notation above, the q-correlation measure jl^ associated with the eigenvalues of a random 
matrix following the distribution F{N, 6, C) is given by the formula: 



'E 



fc.r>l ^-^k.r 



Afe' 



E 

(Afc,r)fc,r>i6L(Af,g) 



n 



N + l -J 

N + 9 



X 



n 

k,r>l 



1 



9 
r\k 



-J 

^k,r 



r\h I ^-^ 



. o-es„ 






(g)Afc , 



.fc>l,r>l 



Remark 3.18. One can check that Proposition 13.121 is a particular case of Proposition 13. 171 
for g = 1. Note that for N,q > 2, the g-correlation measure /i^ is not absolutely continuous 
with respect to Lebesgue measure. This means that the law of the point process of the 
eigenvalues cannot be described with correlation functions, as for determinantal processes. 

In a similar way, we can define the g-correlation measure jl^, associated with the random 
point measure /loo defined by (13. 6p . We obtain the following: 

Proposition 3.19. Let 6 > and C be a probability distribution on C* . Then, the q- 
correlation measure fix, associated with the point measure /ioo defined by (13. 6p is given by: 

^k,r 



H'OD 



E n 



1 
Afc,r! 



9 

r\k 



E- 



TeSa 



.fc>l,r>l 



c^ 

'~'l 



Since we have an explicit expression for the correlations measures, we can expect some limit 
theorems when the dimension A^ goes to infinity. In fact, we have the following proposition: 

Proposition 3.20. Let 9 > Q, q > 1 integer, L a probability measure on C*, (MAr)Ar>i 
a sequence of random matrices such that M^ follows the distribution F{N,9,C). We sup- 
pose that if Z follows the distribution C, then log(|Z|) is integrable, and we define R : = 
exp (E[log(|Z|)]). Then the probability measure j^rj^_^\ ^(n~ +i) I^n (which represents the 
probability distribution of a random sequence of q eigenvalues of Mn, chosen uniformly 
among the N{N — 1)...(A^ — g + 1) possible sequences), converges weakly to the q-th power 
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of the uniform distribution on the circle {z = R}. Moreover, iflog{\Z\) is inW~^^, then for 
< Ri < R < i?2; the restriction of fi^ to the set {{zi,...,Zq) G 0,Vr < q, \zr\ ^ (-Ri,i?2)}, 
converges weakly to the corresponding restriction of jlSo ■ 

Proof. Recall that jr^ is the average of the random measure: 

where {ujj)i<j<n is a sequence containing all the eigenvalues of M^, with multiplicity. Let 
fi, ..., fq be bounded, continuous functions from C to R+. One has: 



I :-- 



which implies, 



for 



fl{Zl)...fq{Zq)dfI^^{Zl,...,Zq) =E 



CI 






iK^jJ 



J-\\fi\U..\\fq\UN''-N{N-l)...{N-q + l)]<I<J, 



J ■=E 



' q / N 

n E/^(^^) 

.r=l \j=l 



or equivalently 



J = E 



f\( f frdf^iM^: 

r-=l V^C 



Now, by Proposition 13.11 and Corollary 



27r 



i- [ frdf^iM^)) -^ (1- I friRe^')dX 

in probability. By applying the bounded, continuous function from M^ to M 

(a;i,...,Xg) (-)■ (|a;i|...|a;g|) A (||/i||oo---||/g||oo), 

one deduces that 

J 1 



m N^^ (27r)5 J [0,2.] 

By using the inequalities (13. lip , one obtains that 
/ 1 



f,{Re''')...fq{Re'^^)dX,...dXq. 



N{N-l)...{N-q + l) N^^ (27r)'. j[o,2.]. 



f,{Re''')...fq{Re'^'')dX,...dX 



(3.1i: 



9' 



which gives the first part of Proposition 13.201 In order to prove the second part, let us 
suppose that < Ri < R < R2 and that / is a bounded, continuous function from C^ to 
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M, vanishing if one of the coordinates has a modulus between Ri and i?2- It is sufficient to 
prove: 



(Afe,r)fc,r>lGL((?) 



X 



n 

fc,T->l 



1 / 9 



Afc,r! \r!A; 



Afc,,- 



ctGEo 






a. 



.k>l,r>l 



O 

'~'i 



jAfc.,. 



N^oo 



0, 



with the notation of the proof of Propositions 13. 15) and [3TT71 By dominated convergence, we 
are done if we check: 



1 / /Q \ Afc 

n 



fc,r>l 



A 



Ms)"] 5/./" 



0". 



.fc>l,r>l 



D"^- 
*-i. 



^Xk,r 



< oo. 



Now, for {Xk,r)k,r>i ^ -^(?)) by the estimates obtained in the proof of Proposition 13.151 

/ \ (i-e)+ 



k.r>l 



< 2qC{e) [1 V sup{A; > 1, 3r > 1, A^,, > 0}]^^"^^^ 



fc,r>l 



where C{6) > depends only on 6. By assumption, if Z is a random variable which follows 
the distribution C, then log(Z) G W^^^, which implies that 

Ck{zeC,\z\^iRi,R2)} 

decreases at least as fast as l/fc'^+^ when k goes to infinity, with a constant depending only 
on C, Ri and i?2- Since / vanishes when a coordinate has modulus in [Ri, R2), one deduces: 



fd 



^k 



»\k,r 



<c n 

fc>l,r>l 



A;(A;-l)...(A;-r + l) 



.fc>l,r>l 

where C can depend on £, i?i, R2, f, q. Hence, one only needs to check that: 



E n 



1 fee, 



XkJ V k^ 



< 00, 



where /3 > g + 1 and where by convention, the binomial coefficient C^. is equal to zero for 
k < r. Now, for (Afc,,.)fc,r>i G L{q)- 



and 



C^< fc"^ 
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Hence it is sufficient to have 

(Afc,.)fc,.>ieZ(q)fc,r>l 

Since 

k,r>l 

one necessarily has Xk,r = for r > q. Hence, in any case, 
for 7 = g — /3<— 1, and Proposition 13.201 is proved if 

j2 n '^^^'■" < ^• 

{^k,r)k.r>i(^L{q)k,r>l 

This last estimate is easy to check, by replacing L{q) by the (larger) set of families {Xk,r)k,r>i 
of nonnegative integers such that A^^^ = for r > q and Xk,r < Q' in any case. D 

4. The point process of eigenangles and its scaling limit in the unitary case 

In this section we do a more precise study of the point process of eigenvalues, in the 
unitary case. In particular, we obtain a scaling limit for the eigenangles if they are properly 
renormalized. This limit is dominated by the large cycles of permutations, hence, we can 
expect almost sure convergence results if we consider virtual permutations. More precisely, 
let {Misf)i\f>i be a sequence of random matrices following the distribution P(cxd, 9, £), where 
6 > and £ is a probability distribution on the unit circle. Recall that the measure ^{M^), 
representing the point process of eigenvalues can be written as: 



^[M^) = Y^ i,^^^^>o Y, ^^ 



"^=1 a;'iV.,n = v ^^ 



Here lN,m. is the cardinality of the intersection CiVm of C^ and {1, ..., A^}, and 



ZN,m — _[_[ 



Zj, 



where (Cm)m>i is the partition of N* (ordered by increasing smallest element) given by the 
cycle structure of a virtual permutation following the Ewens(6') measure, and {zj)j>i is an 
independent sequence of i.i.d. variables with law C. Since all the eigenvalues are on the 
unit circle, it can be more practical to consider the eigenangles. The corresponding point 
process can be described as follows: one takes the point x if and only if e*^ is an eigenvalue 
of Mat, with a multiplicity equal to the multiplicity of the eigenvalue. Note that this process 
is 27r-periodic, and the corresponding random measure t{Mn) can be written as: 

oo 

t{M^) = Y lu,„.>o Yl '^- 

m=l x=ejv,m(mod. 27r/i]v,,„) 
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where 

©TV.m = -J J^ ^Og{Zj) 

is real and well-defined modulo 2ti /lN,m, for l^^rn > 0. Note that the construction of the 
measure t{Mi\i) implies immediately that t{Mn){[0,2it)) = N, in other words, the average 
spacing of two consecutive points of the corresponding point process is equal to 271 /N. If we 
want to expect a convergence for N going to infinity, we need to rescale t{Mn) in order to 
have a constant average spacing, say, one. That is why we introduce the rescaled measure 
T]\f{Mj^), defined as the image of t{Mn) by the multiplication by N/2tt: the corresponding 
point process contains the points x such that e^'^'^^l^ is an eigenvalue of Mat. One checks 
that 

oo 

TAr(M7v) = 2_^ lyiv,„>o2^'^(7iv,m+fc)/yiv,m' 
m=\ fceZ 

where y^^^ra '■= ^N,m./N is the m-th renormalized cycle length of the permutation ctat associated 
with Mat, and: 



^^'--■=^ E 1°^ 



2.vr . ^ -"-^"- 

is well-defined modulo 1. Now the general results on virtual permutations (see [H] or Section 
[2]) imply that for all m > 1, yN,m converges a.s. to a random variable ym, where {ym)m>i is a 
GEM process of parameter 9. Then, if •yN^m is supposed to be equal to zero, one can expect 
that the measure t^{M^) converges to 

oo 

'Too{{Mn)n>i) ■■= Yl ^h/y^- 
m=l feeZ 

Of course, one needs to be careful since Too((M7v)Ar>i) has an infinite Dirac mass at zero. 
Moreover, the condition ■yN,m = is satisfied if and only if Cis 6i, the Dirac measure at one. 
That is why we state the following result: 

Proposition 4.1. Let (MAr)iv>i be a sequence of random matrices, which follows the distri- 
bution F{oo,6,6i), for 6 > 0. Then, with the notation above, the random measure tn{Mn) 
converges a.s. to the random measure Too{{M n) n>i) , in the following sense: for all contin- 
uous functions f from M to M+, with compact support, 

fdTN{MN) -^ /"/rfroo((M^)7v>i) 

almost surely. 

Remark 4.2. In Proposition 14. ![ the matrix Mn is simply a permutation matrix, associated 
to a permutation a^ G ^n- Moreover, the positivity of / is needed in order to deal with the 
infinite Dirac mass at zero of Too((M7v)Ar>i). 

Proof. If /(O) > 0, the integral of / with respect to tj^{Mj^) is greater than or equal to /(O) 
times the number of m such that yN,m > 0, i.e. the number of cycles of a^. Since this 
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number tends a.s. to infinity with A^, we are done. Then, we can suppose /(O) = 0, which 
imphes: 



p oo 



m=l fceZ\{0} 

and 

» CO 

"^^ m=\ fceZ\{0} 

where we take the convention f{k/yN,m) = for yN,m = 0. Now, if the support of / 
is included in [—A, A] and if Sm is the supremum of yN,m for all A^ > 1, it is clear that 
f{f^/yN,m) = f{k/ym) = for all k,m such that \k\ > Asm- One deduces that 






and 



/"/droo((M^)iv>i)= Yl /(^/^-) 



{m,k)(^S 

where S is the set (independent of A^) of couples of integers [m, k) such that m > 1 and 
< |fc| < Asjn- Now 

Ysm<00 
m>l 

almost surely (see (I3.3P ). and then S is a.s. a finite set. Since fik/yN^m) tends a.s. to 
f{k/ym) for all k,m (recall that / is continuous), we are done. D 

The a.s. weak convergence given in Proposition 14.11 cannot be directly generalized if the 
law £ is not the Dirac measure at one, because ■yN,m changes in a non-negligible way for 
each N G Cm- However, one can expect a weaker convergence. A natural candidate for the 
corresponding limit distribution would be the law of a random measure defined by: 



X]X^'^(fc+X.n)/?/m' (4-1) 



m=l 



where we recall that {ym)m>i is a GEM(6') process, and where (Xm)m>i is an independent 
i.i.d. sequence in M/Z. The distribution of Xi needs to be, in a sense which has to be made 
precise, close to the distribution of: 

1 ^ 

where X is a large integer (recall that {zj)j>i is an i.i.d. sequence of random variables with 
law £). From the following result, we deduce a good candidate: 

Lemma 4.3. Let C be a probability measure on U, and for k > 1, let C*^ be the multiplicative 
convolution of k copies of C Moreover, let r{C) be the infimum of the integers r such that 
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L is carried by the set Ur of r-th roots of unity. Then for all sequences {dk)k>i of strictly 
positive integers, tending to infinity with k, the probability measure: 



fc+dfe-l 
d~ 



1 ^ ^*. 



p=k 

converges weakly to C* , where C* is the uniform distribution on Uj. if r{C) < oo, and the 
uniform distribution on U ifr{C) = oo. 

Proof. Let us define, for all g G Z, and for all probability measures Ai on U: 

M(q) ■= f z'^dM. 

Jv 

Moreover, let us set, for A; > 1: 

fe+dfc-l 
p=k 



Then, one has: 

k+dk-l 
p=k 

which is equal to one if C{q) = 1 and which tends to zero when k — )■ oo, if C{q) ^ 1. Now, 
C{q) = 1 if and only if q is divisible by r{C), for r{C) < oo, and if and only if g = 0, for 
r(£) = oo. Hence: 

^*(?) = l£(g)=0' 

and finally, 

Mk{q) -^ C*iq). 

D 

Because of Lemma [4.31 the law D{C) of xi, in equation (14. ip . is chosen as follows: if r{C) < 
oo, it is the uniform distribution on the classes {0, 1/r, ..., (r — l)/r} modulo 1, and if 
r(£) = oo, it is the uniform measure on M/Z. 

Remark 4.4. If r{C) is finite, it is the smallest integer r > 1 such that M^ G T-Lri^) a.s., for 
Mn following the distribution P(iV, 6, C). 

We are now able to state our convergence result for any distribution C on U: 

Proposition 4.5. Let (Mjv)jv>i be a sequence of random matrices, such that for all N >\, 
Mjy follows the distribution F{N, 9, C), where 6 > Q and where C is a probability measure on 
U. We suppose that C satisfies one of the two following conditions: 

• the measure is carried by U^ for some integer r > 1 . 

• there exists f > 1 such that for e > small enough, and for all arcs A in U of size 
e,C{A)<\\og{e)\-r 
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Let Too{0,jC.) be the random measure defined by: 

oo 
m=l keZ 

where {xm)m>i is a Poisson-Dirichlet process of parameter 6 , and {Xm)m>i is an independent 
sequence of i.i.d. random variables on M/Z, following the distribution D{C) defined above. 
Then, with the previous notation, the distribution of the random measure tn{Mn) converges 
to the distribution of the random measure Too{6, C), in the following sense: for all continuous 
functions f from R to M+, with compact support, 

• if f{^) > and r(£) < oo, then 



a.s., and for all A > 0: 



P 



/ fdT^{e,C) = oo 

Jm 



fdTN{MN)<A 



when N -^ oo; 

if f{^) = or r{C) = oo, then: 



[ fdT^{e,c)<oo 



a.s., and 



in distribution. 



fdTN{MN) 



N^oo 



fdT^{9,C) 



Remark 4.6. In Proposition 14. 5[ the Poisson-Dirichlet distribution can be replaced by a GEM 
distribution, since it does not change the law of the random measure t^{9, C). 

Proof. Without changing the laws of the random measures tj^^M^) and t^{0,C), one can 
suppose that {M]si)n>i follows the distribution P(oo, 6, C) and that, with the notation above: 



m=l keZ 

where (Xm)m>i is supposed to be independent of {ym)m>i- Recall that 



m=l kez 

where, conditionally on {yN,m)m>i, (7Af,m)m>i is a sequence of independent random variables 
on M/Z, and for all ?7i > 1, 'jN,m has the same law as -^f;;-, where Z follows the distribution 
£*NyN,n, = £*iN,,n_ Lg^ ^g j^q^ supposc /(O) > and r(£) < oo. Since Xm is a.s. equal to 
zero for infinitely many m (because D{C){{0}) > 0), Too{6,C) has an infinite Dirac mass at 
zero, and then a.s.: 

/ fdT^{e,C) = oo. 

Jr 
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Moreover 

oo 

fdr^iMr,) > /(O) 5^ l,^,„>o,7^,„=o. 

m=l 

Now, if (^p,g)p,g>i is a family of independent random variables, independent of a^ (the 
permutation associated with M^), and such that for all p,q > 1, Zpg has distribution C*p, 
then 

oo oo <^N,p 

/ , lj/]V,m>0,7iV,m=0 = 2_^ 2-^ 1^P,? = 1 

m=l p=l (jr=l 

in distribution, where ajsi^p denotes the number of m such that /jv,m = P (i-e. the number of 
p- cycles in (Tn)- Therefore: 

„ oo cjv,p 

•^^ p=l g=l 

where ^ denotes the stochastic domination, and where {cN,p)p>i, defined by the Feller cou- 
pling as in (13.71) . is supposed to be independent of (^p,g)p,g>i- Since for all p > 1, cn,p 
increases a.s. to 6p, Proposition 14.51 is proved for /(O) > and r(£) < oo, if we show that 
a.s. , 

oo bp 
p=l q=l 

and a fortiori, if we prove that 

oo 

/ , '^3q<bp,Zp,q=l = OO. 
p=l 

Since all the variables (&p)p>i and (^p,5)p,5>i are independent, by the Borel-Cantelli lemma 
it is sufficient to have 

oo 

J]P[3g<6j„Zp,, = l] = oo. 
p=i 
Now, if for p > 1, we define 

p{p) = r^({i}), 

we have: 

P[3g < bp, Zp,, = l\bp] = l-[l- P{p)]'% 
and, since bp is a Poisson random variable of parameter 6/p, 

P[3g<6p,Zp,, = l] = l-e-^^(P)/^ 
and then we are done if we prove: 

5^P(p)/p = oo. (4.2) 

p>i 

Now, by Lemma [4.31 
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converges weakly to the uniforin distribution on Ur{c), for k — )■ oo. Since all the measures 
involved here are supported by the finite set Vr(c), one deduces that 



2k+i_i 



i E P(P) .-1 Vr(£) 



=2fe 



fc— >oo 



which implies fl4.2p . We can now suppose /(O) = or r{C) = 00. In order to prove 
Proposition 14.51 in this case, we need the following result: 

Lemma 4.7. Let 6 > 0, r & W , and let {lN,m)m>i be the sequence of cycle lengths (ordered 
by increasing smallest elements, and completed by zeros) of a random permutation in T,^, 
following the Ewens measure with parameter 6. Then, there exists a function \l/r,6» from W 
to R"*", tending to zero at infinity, and such that for all strictly positive integers {lm)i<m<r, 
(Oi<m<r, N, satisfying: 

m=l 
m=l 

Vm G {1, ...,r}, 3g G N*, q^ < /„, C < (g + l)^ 
m G {1, ...,r}, /^,„ = ^J ^ (1 _ ^^^^(^)^ 1 + ^^^^(^)) . 



one has 

P[\ 

PJVmG{l,...,r},/^,^ = /J 
Proof. For < ?72 < r, let us define: 



/ Jp-' 
p=i 



^m — 2^ ^p' 



p=l 

By using the Feller coupling, one obtains the following expression: 

/m G {1, ..., r}, lN,m = ^m]=Yl 



L, 



9 



n 



N -p 



N-Lm + e ^i-i. N-p + e 
1+1 



n 



iV-p 



P=Lm- 
T 

n 







A^-L, 



Hence, if L^ > L'^ 



m G {1, ...,r},/ 



Ar,m 



/' 1 




N -p 

N-p + e 



TT AT-L'^ 



\m=l 



N-Lr. 



For 1 < m < r, there exists an integer g > 1 such that q"^ < Im, I'm < (q + ^Y ^ ^-iid q < vN, 
since ImJ'm < ^- Hence, \l.m — I'^l < 2g < 2\/N, and \Lm — L'^\ < 2r\fN. Moreover, 
A^ — Lm, A^ — L'^ > N'^l^ by assumption. Therefore, one has the majorization: 



P[Vme {l,...,r},/ 



N,m 



°^VP[VmG{l,...,r},/^,, 
This imphes Lemma 14.71 for: 



< 2r\/N log 1 + 



Ar2/3 



+ r log 1 + 



2rViV 



2rVN 



N^/sJ 



2rVN 
iV2/3 



+ 



TV 



1. 



D 



A consequence of Lemma 14.71 is the following: 



Lemma 4.8. Let 6 > and let C be a probability law on the unit circle. Let {lN,m)m>i 
be the sequence of cycle lengths (ordered by increasing smallest elements, and completed by 
zeros) of a random permutation in Sjv, following Ewens measure with parameter 6' > 0. Let 
{,lN,m)m>i be a sequence of random variables on M/Z such that conditionally on {lN,m)m>i, 
i'yN,m)m>i are independent and for all m > 1, '')N,m has the same law as -^^, where Z 
follows the multiplicative convolution C*^N,m of l^^rn copies of C Then, for all r > 1, the 
family {lN,m/N,'yN,m)m<r converges in distribution to {ym,Xm)m<r, where {ym)m>i is a GEM 
process of parameter 6 , and {Xm)m>i is an independent sequence of i.i.d. random variables, 
with law D{C), defined above. 

Proof. Let qN,m be the integer part of the square root of lN,m- For all continuous functions 
g from (M/Z)'' to M+, one has: 



ngih 



)l(9. 



N,in)m<r ) \qN,m)m< 



:r]= / g{\ogZi/2i7r,...,\ogZr/2iTT)dAf(^ 

Jw 



(giV,m)m<rV^l' ' ' ' ' ^''J' 



where A/(gjvm)m<r is defined as the conditional expectation of a random measure, more pre- 
cisely: 



M 



{QN,mjm<T 



E 



C 



''jV.i 



m=l 



(q 



N,m)m<r 



Let us suppose that qN,m > N^^^ for all m < r and: 

r 



m=l 



By Lemma [4. 7[ the conditional law of {lN,m)m<r given {qN,m)m<r is a linear combination of 
Dirac measures on the r-uples {lm)m<r of integers such that qjf m — ^m < {qN,m + ^Y, and the 
quotient between two coefficients of this combination lies on the interval (1 — \l/r,e(A^), 1 + 
\l/r,6i(A^))- One deduces that 



(1 - ^rAN)) (8)C < ^i,..mU^r < (1 + ^rAN)) (g)/:; 



qN,n 



m=l 



m=l 
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where, for g > 1: 



q^+2q 



CL 



2q + 



1 -^ ' -i 

-^ H ^*' 

+ 1 ^^ 



Hence: 



noiir 



)l(?iV,r 



N,m)'m<r) \hN ,m j ni<r 



C [ gi\ogz,/2tn,...,\ogZr/2tn) HC'^^Jdz 



for 



Now, by Lemma H73l Cg converges to C* when q goes to infinity, hence, there exists a function 
Kg from M_|_ to R_|_, tending to zero at infinity, and such that for all q > qo ^ ^, Q integer: 



/ g{^ogzi/2iTi, . . .,\ogZr/2i7r) T\ C'g^^{dZm) 



- / g{\og Zi/2i'K,..., log Zr 12171) \\C*{d. 
Therefore, if N is large enough, if qjss^m > N^^'^ for all m < r and if 

then: 



< ^^(^o). 



m=l 



E[^((7. 



N,m)m<j' 



')\iQN., 



m)m<r 



e[A,,B, 



for 



Ai = (l-^,,e(iV)) 



and 



5i = (1 + ^rAN)) 



-Kg{N^I'')+ [ g{logzj2tn,...,logZr/2t7i) Jl D{C){d2 

J^^ ra=l 

Kg{N^'^)+ [ g{l0gZi/2t7l,...,\0gZr/2l7r) U D{C){dZrr. 



Now, let / be a continuous function from [0, 1]^ to M+. One has: 

W{{(il,JN)^<r)g{{lN,^)m<r)\ G [^2,^2], 

for 

A^ = A,E[f{iqlJNU^r)te] 
and 

B2 = B,E[fi{ql,jNU<r)le] + ||/||oo||^||ooP[^^ 
where 



S := {Vm < r, qN,m > N'/\ N - J2i^N,m + if > iV'^'} 



m=l 



and S'^ is the complement of the event S. If we suppose that {lN,m)N,m>i are the cycle lengths 
associated with a virtual permutation following the Ewens(^) measure, then for all m < r, 
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q%^/N tends a.s. to ?/„, for a GEM(^) processes (ym)m>i, and the event E'^ holds for finitely 
many values of N . Therefore: 



N^oo 



and 



One deduces that 



P[£1 -^ 0. 

Af-5-oo 



Hf{{(&jN)ra<r)g{{lN,m)m<r)] -~^ ^f {{V m) m<r) 9 {{Xm) m<r)] 

(recall that {ym)m<j- and {Xm)m<r are supposed to be independent). Now, 



lN,m V 



N'tri 



N N 



< 



N 



hence, by the uniform continuity of / and the boundedness of g: 

Hf{lN,m/N)m<r)g{{lN,m)m<r)] > ^[f i{ym)m<r)g{{Xm)m<r)], 

which proves Lemma [4.81 
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Now, let / be a continuous, nonnegative function with compact support, and $ a contin- 
uous function from M+ to [0, 1], such that $(a;) = for x < 1 and $(x) = 1 for x > 2. Let 
us suppose /(O) = or r{C) = oo. We first remark that for all integers r, s > 1: 



m=l fceZ m=l fceZ 

in distribution. This is a consequence of Lemma 14.81 and the fact that 

^sy)J2m + x)/y] 

fcez 

is continuous with respect to {y,x) ^ [0, 1] x (M/Z). Now, if / is supported by [—A, A] (for 
y4 > 0), one has the inequalities: 



and 



J] /[(A; + lN,m)/yN,m] < {2A + 1 
J2flik + Xrn)/ym]<{2A+l 



kez 



Therefore, for A G M+: 

-^ g«^Em=llj/iV,m>''Efc6z/[('^+')'JV.™)/S'JV,m] 
< ]g p«^Em=l*(«S/JV,m)Efc6z/[('=+7iV,7n)/yiV,m] 



E 

-E 



= «-^ Em=l EfcsZ /[C^ + Xm)/!/™] 



o*-^ Em=l *(«2/™) EfcgZ /[C^ + Xm)/!/™] 



+ (2A + 1)A||/|U Yl (^[0 < ^^-- ^ 2/s] + P[|/^ < 2/s]) . 



m=l 
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The first term in the right-hand side of this inequality tends to zero when A^ goes to infinity. 
Hence: 



hm sup 

Af-s>oo 



E 



«^Em = llaiVm>0Efc6z/[(^+7iV,m)/l/iV,m] 



-E 



3«^ Em=l EfceZ /[(fc+Xm)/i/m] 



< 



m=l 



{2A + l)A||/||oo J2 (limsupP[i/;v,n^ < 2/s] + P[y„ < 2/s] 



Af-s>oo 



Now if we assume that {yN.m)N,m>i is the family of renormalized cycle lengths of a virtual 
permutation following the Ewens(6') measure, we obtain, by Fatou's lemma: 



liminfPfi/iVm > 2/s] >E 
Since yN,m tends to ym for N going to infinity: 



lim inf 1 



yN,m>'i/s 



limsupP[y^,^ < 2/s\ < P[y™ < 2/s]. 



By taking s — )■ cxd, we deduce that: 



E 



o*-^Em = llyjV,m>"Efcez/[(^+TJV,m)/yiV,m] 



-E 



=«->* Em=i Efcez /[(*:+Xm)/yn 



(4.3) 



tends to zero when A^ goes to infinity. Now, for all x ^ [O? 1) ^"^^ y ^ (0) !]• 



Y,fU + x)/y]<{'iA + i)\\f\UiAy>i. 



and if /(O) = 0, X = 0: 



Therefore, if 7Af,m is identified with the unique element of [0, 1) in its congruence class modulo 
1, then: 



E 



Y^ ^yN,n.>oYf[i.k + lN,m)/y_ 



'N,m\ 



m=r+l 



<(2A+1)||/|U 5^ [P(y7V,™ > 0, < 7jv,,n < A2/iV,n^) 
m=r+l 

+^{yN,m > 0, 1 - 7Ar,m < AyN,m) + '^{lN,m = 0, AyN,m > 1)] 



(4.4) 



if /(O) = 0. In fact, this inequality remains true for /(O) > 0. Indeed, in this case, one has, 
by assumption, r(£) = oo, and then by the conditions given in Proposition 14. 5[ C has no 
atom. Since 



E 



Yl lyiv,™>oXl.^[(^ + ^^'"^)/^^'' 



,m=r+l 



oo 



< (2.4 + 1)1 1 /I loo Y \^(yN,m > 0,7JV,m < AyN,m) + HyN,n. > 0, 1 - 7^,^ < AyN,m)] 



m=r+l 
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we also have f l4.4p . Now, for all p G (0, 1): 

P[l/^,^ > 0, < ^N,m < AyN,m] = P[Z/7V,n^ > 0, < ^N,m < P™] + P[l/7V,m > p'^/A]. 

Since conditionally on {yN,m > 0}, the law of 'jN,m is a convex combination of iterated 
convolutions of £, it satisfies the same assumptions as £, with the same underlying constants. 
One deduces that there exist f > 1, rriQ > 1, C > 0, independent of A^, such that for m > niQ-. 



P[2/iv,m > 0, < 7^,^ < p"^] < Cm-\ 



(4.5) 



Now, if the normalized lengths of cycles {yN,m)N,m>i are associated with a virtual permuta- 
tion following the Ewens(6') measure, and if for all m > 1, Sm is the supremum of yN,m for all 
A^ > 1, then the expectation of Sm decreases exponentially with m (this result is contained in 
our proof of (13.41) above). One deduces that there exist p G (0, 1) and rrii > 1, independent 
of A^ and such that for m > ttt-q: 



nyN,m > P'^/A] < p™. 



From (US]) and (gSD, 



(4.6) 



sup V" f'[yN,m > 0, < -fN,m < AyN,m\ > 0. 

N>1 ^-^1 r^°° 

— m=r+l 

In a similar way, one can prove: 



sup y~] F[yN,m > 0, 1 - 'jN,m < Ay 

N>1 ^-^., 
— m=r+l 



N,m\ 







and 



Therefore, by (14.41) : 



sup y~] P[7Ar,m 

N>1 ^-^., 
— m=r+l 



0, AyN,m > 1] ^ 0. 



supE 

N>1 



Yl '^yN,m>oYf[{k + -fN,m)/yN,r^ 



_m=r+l kGl 

By replacing 'yN,m by Xm and yN,m by ym, one obtains: 



0. 



E 



Yl i?/->oX]'^[(^+^™)/^" 



m=r+l 



(4.7) 



(4.8) 



since the law of Xm satisfies the same assumptions as £ and the expectation of ym decreases 
exponentially with m. Since the quantity given by (14.31) tends to zero when A^ goes to 
infinity, one easily deduces from (14. 7p and (14.81) that it is also the case for: 



E 



«^Em = llaiV„>0Efcez/[('=+7iV,m)/yiV,m] 



E 



3«^ Em=l EfceZ /[(fc+Xm)/j/m] 



which proves Proposition 14. 5[ 



D 
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5. The uniform case 

5.1. Eigenvalues distributions and correlation measures. In this section, we focus on 
the uniform case, i.e. the case where for A^ > 1, the random matrix M^r follows the law 
P(A^, 6*, £), where 9 > and C is the uniform distribution on the unit circle. Let {yN,m)N,m>i 
be the family of renormalized cycle lengths corresponding to a virtual permutation following 
the Ewens(^) measure, i/m the limit oiyN^rn for N going to infinity, and (Xm)m>i a sequence of 
i.i.d. uniform random variables on [0, 1), independent of {yN,m)N,m>i- The random measure 
tn{Mn) has the same law as 



and Too (6*, £) is equal to 



m=l feeZ 



Too :- 2_^ 7 ,^(Xm+k)/y„^ 
m=l kel 

in distribution. This description of the law of Ti<^{M]\f) and Too{6,C) implies the following 
remarkable property: 

Proposition 5.1. Under the assumptions given above, the distributions of the random mea- 
sures tn{Mn) and Too{6,jC.) are invariant by translation. 

Proof. It is enough to prove it for tn and Too. The image of tn by a translation of A G M is 
equal to 

oo 

/ . ls/iV,m>0 / .^(Xm+fc)/yiV,m) 

m=l kez 

where x'm is the fractional part of Xm + AyN,m- Now one easily sees that conditionally 
on {yN,m)m>i, the sequence {x'm)m>i is an i.i.d. sequence of uniform variables on [0, 1), as 
(Xm)m>i- This implies the invariance by translation of f^. For foo, the proof is exactly 
similar. D 

The main interest of the introduction of the measures tn and Too is the following: we cannot 
expect an a.s. convergence of the random measure tm{Mm) when N goes to infinity, even 
if {Mis[)n>i follows the distribution P(oo, 6*, C) because, with the notation of the previous 
section, the "shifts" '-)N,m (which are equal to Xm in distribution), do not converge a.s. when 
N goes to infinity; however, since we take the same variables Xm in the definitions of f^v and 
Too, we can expect such a convergence for these new measures. More precisely, we have the 
following: 

Proposition 5.2. Almost surely, with the notation above, the random measure f^r converges 
vaguely to f^o, which is locally finite. 

Proof The fact that foo is locally finite is a consequence of Proposition 14.51 (one has r(£) = 
oo). Now, let / be a continuous function from M to M, with support included in [—A, A] (for 
some A > 0). One has: 



« oo 

/ fdfN = J2Yl lj/iv,m>o/((Xm + k)/yN,, 

•'R rr, — ^ l.r-'TI 
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and 






Let Sm be the supremum of yN,m ior N > 1. For \k\ > A + 1 or Asm < Xm < 1 — Asm, one 
has: 



l?/iV,m>o/((Xm + k)/yN,m) = f {{Xrn + k)/y^) = 0, 



and since 



E 



m>l 



< CXD, 



there exists a.s. a (random) finite subset S* of N* x Z such that: 

fdfN= Yl '^yN,^>of{iXm + k)/yN,n 

^^ (m,k)eS 

and 

fdfoo= V f{{Xm + k)/y„,). 

Since a.s., for all m > 1, fc G Z: 



(m,fc)e5 



lw™>o/((Xm + fc)/yiV,m) ^7-> /((Xm + ^)/ym), 



one deduces Proposition 15. 2[ 



D 



One can now study the correlation measures associated with the point processes studied 
above. More precisely, if 

where {xj)j>i is a sequence of real numbers, one can define, for all g > 1 the random measure 
on W: 

rf(M^)= E ^^n^-,^..)- 

Similarly, one can define TcS (6',£), r^ and r^ . Of course, one has the following equalities 
in distribution: 



r].^'(M^) = r 



\i\ 



N 



and 



r!^{0,C) = f^^. 



Proposition 15.21 can be generalized as follows: 



J-?] 



M 



Proposition 5.3. For all q > 1 the random measure f™ converges vaguely towards f^ , 
which is locally finite. 
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Proof. The fact that f^ is locally finite is a consequence of the local finiteness of foo- Now, 
let / be a continuous function from M^ to M, with support included in [—A, A]^ (for some 
A> 0). One has: 

/ fd^N = Yl lw,™i,...,2/iV,™,>0/((Xmi + ki)/yN,rm, •••, {Xm, + kq)/yN,m,) 

and 

/ fdf^^^ = 53 /((Xmi + kl)/yrm, •-, (Xm, + ^g)/Z/mJ 

Let Sm be the supremum of yN,m ioi N > 1. If for some j < q, \kj\ > A+1 or Asm^ < Xmj < 
las: 

lyiV.™i,...,J/iV.m,>o/((Xmi + ki)/yN,rm, •-, (Xm, + kq)/yN,m,) = 



1 — Asmj , one has: 



and 

fiiXrm + ki)/ym„ ..., (Xm, + kg)/ymg) = 0. 
Since 



< oo, 



E Y.^ 

m>l 

there exists a.s. a (random) finite subset S* of N* x Z such that: 

/ /c^^Jv^ = Yl lyiV,™i,...,J/iV,™,>0/((Xmi + ki)/yN,rm, •-, (Xm, + kg)/yN,m,) 

and 

/ /C?^i^^ = X] fiiXrm + ki)/yrm, •••, (Xm, + kq)/ym,) 

Since a.s., for all mi, ..., rrig > 1, fci, ..., fcg G Z: 

l|/iV,™i,...,to,^,>0/((Xmi + ki)/yN,rm, •-, (Xm, + kq)/yN,m,) 

tends to 

/((Xmi + ki)/ym^,..., {Xm, + fc<?)/l/mj 

when A^ goes to infinity, one deduces Proposition 15.31 D 

Now, for q > 1, we define the "g-correlation measure" fjy associated with the point process 
T]\f{M]\i) (or equivalently, tat), as the average of the measure r]^'(Mjv) (or f^^)- Similarly, 
we define the measure fs} as the average of t^{6,C) of f^ . Both r^' and f^ are positive 
measures on W. From the convergence of fjy toward fti^ , one deduces the following result: 

Proposition 5.4. For all q > 1, the measures fjy (N > I) and f^ are locally finite, and 
tJy converges vaguely towards f^ when N goes to infinity. 
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Proof. Let / be a continuous function from M to M, with compact support. One has, 



fdf^^^ 



< 



E 



1. 



°° / ■ "{\kj\<A+l)j<g,{xm^t(Asn.j,l-Asmj)) 

(mi,ki)^...^(mq,kq)eN*xZ 



■3'lj<q 



which is independent of A^ and has a finite expectation. Indeed, this quantity is bounded by 
a constant times 

/ , (lxm^{As„,l-As„)j 
m=l 



and 



E 



E(i 



Xtti^ (-^^771)1 -^^ra) } 



\ra=^\ 



< E n{Xrn,^iAs„l~As,))^^^] 



mi,...,mq'>l 



ni>l max{mj,l<j<(jr}=m 

< 2A ^ m«E[s„] < oo. 



m>l 



since the expectation of Sm decreases exponentially with m (see the proof of fl3.4p ). Almost 



surely, by proposition 15. 2[ 



fdfl^^ 



N^oo 



and then one obtains, by taking expectation and applying dominated convergence: 

fdf^N ;r^ [ fdf!^\ 
where all these integrals are finite. This proves Proposition 15.41 



D 



For q = 1, the g-correlation measure of a subset of M is simply the average number of points 
lying on this set. It can be very simply expressed: 



^[1] 



^[1] 



Proposition 5.5. The 1- correlation measures Tj^ (N > 1) and t^ are equal to Lebesgue 
measure on M. 

Proof. Let / be a nonnegative, continuous function from M to M. One has: 



E 



fdz 



N 



iVN, 



in)in>l 



^yN,m / f{z)dz 
m>l -^^ 

f{z)dz, 



^w 



which proves Proposition 15.51 for r|r . The proof for f^ is similar. 
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We remark that for g > 2, the correlation measure fj^ is not absolutely continuous with 
respect to the Lebesgue measure. Indeed, for all integers I, 1 < I < N, the probability that 
the point process associated with fjv has two points separated by an interval of exactly N/l is 
not equal to zero (this event holds if the corresponding permutation has a cycle of length /). 
Similarly, for g > 3, f^ is not absolutely continuous with respect to the Lebesgue measure, 
since the point process associated with foo has almost surely three distinct points a;, y, 2; G M 
such that y — X = z ~ y. 

~\2] 

However, despite the fact that rj^ is not absolutely continuous with respect to the Lebesgue 

~[2l 

measure for all iV > 1, its limit rio is absolutely continuous with respect to the Lebesgue 
measure. More precisely one has the following: 

Proposition 5.6. The measure fU on R^ (which depends on the parameter 9) has a density 
p with respect to the Lebesgue measure, which is called "2- correlation function", and which 
is given by: 

p{x,y) = (j)0{x-y), 
where the function (j)^ from M. toM. is defined by: 



9+1 



x^ 



aeN*,a<\x\ ^ 

Proof. One can write: 

m^m'>lfc,fc'eZ 

+ 2^ 2^ '^((Xm+fc)/j/m,(Xm + fc')/S/m)5 

m>lfc^fc'eZ 

which implies that for all nonnegative and continuous functions / from M^ to 

-1 /"I 



E 



fdft^ 



(y. 



m)m>l 



= ^ ^ I I f{{k + x)/ym,{k' + x')/ym')dxdx' 
+ ^ ^ / f {{k + x)/ym, {k' + x)/ym) dx 






m>laGZ\{0} 

The expectation of the first term of (15.11) is equal to 



f{x,x + a/yrr,)dx. 



(5.1: 



1-E 



E 

.TO>1 






f{x,x') dxdx' 



9 



9 + 1 



f{x,x')dxdx'. 



In order to compute the expectation of the second term, let us fix m and a, and let us 
denote by dm the density of the probability distribution of ym (recall that dm{x) = for all 
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X ^ [0, 1]). One has: 



¥.[y^f{x,x + a/ym)\= / tdm{t)f{x,x + a/t)dt 

Jo 

.2 



a 



\u\ 



and then: 



E 



Vm / fix,x + a/ym)dx 






:dm{ci/u)f{x, X + u)du 



:drn{a/{x' — x))f{x,x')dxdx'. 



Finally, the expectation of the second term of (15.11) is: 



aeZ\{0} rn>\ ' 

which proves Proposition 15.61 with 



\x — x\ 



-dm{o,/{x' — x))f{x, x')dxdx' 



9 



Y+ u:pZl^'Zl^™^^/l^l)• 



a>l m>l 

Now, for all continuous functions g from [0, 1] to M+, one has: 



E 



XI ymdiy^ 



m>l 



ngivAi)], 



where conditionally on {ym)m>i, the random index M is chosen in a size-biased way, i.e. 
M = m with probability y^- By classical properties of GEM and Poisson-Dirichlet processes, 
yM is equal to yi in distribution, and then its density at x G [0, 1] is 9{1 — x)^~^. Hence, one 
deduces that 

y dm{x) = 9 (1 — x^^^gix) dx, 

m>l J^ 



xg[x) 



and then, for almost every x E [0, 1], 

m>l 

which implies Proposition 15.61 



X 



9{i-xy-^ 



X 
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5.2. The smallest eigenangle. Another interesting problem about the point process asso- 
ciated with tn{Mm) is the estimation of its smallest positive point. This point corresponds 
(after scaling the eigenangle by A^) to the first eigenvalue of M]si obtained by starting from 
1 and by turning counterclockwise on the unit circle. Its distribution has a limit by the 
following result: 

Proposition 5.7. With the notation above, the smallest positive point corresponding to the 
random measure f^ tends a.s. to the smallest positive point corresponding to foo- 
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Proof. One has, for all m > 1, 

if by convention, Xm/yN,m = +00 for yN,m = 0. One deduces that for all mo > 1, 

lim sup mi{xm/yN,m, 1 < 
mi{xm/ym, 1 < m < mo} 



Xm/yN,m ~ ' Xm/yn 



lim sup mi{xm/yN,m, m > 1} < lim sup mi{xm/yN,m, 1 < m < mo} 



and then, by taking mo -^ 00: 

lim sup mi{xm/yN,m,m > 1} < mi{xm/ym,m > 1}. 

On the other hand, for all A > 0, there exists a.s. m-i > 1 such that Xm/sm > A for all 
m > mi, which implies that Xm/yN,m > A and Xm/ym > ^- Consequently 

lim inf mi{xm/yN,m,m > 1} > A A lim inf mi{xm/yN,m, 1 < m < mj 

> A A mi{xm/ym, 1 < m < mj 

> A A inf {xm/2/m, "^ > !}• 
By taking A — )■ 00, one obtains: 



and finally. 



lim inf mi{xm/yN,m,m > 1} > mi{xm/ym,m > 1}, 



mi{xm/yN,m,rn > 1} — > mi{xm/ym,m > 1}, 



which proves Proposition 15.71 D 

One immediately deduces the following: 

Corollary 5.8. The smallest positive point of the random measure tj^{M]\}) converges in 
distribution to the smallest positive point of Tao{0,C). 

At this stage, one naturally seeks for the explicit computation of the distribution of the 
smallest positive point of Too{6,C). We can remark the similarity between this process and 
the limit point process (a determinantal process with sine kernel) obtained from the scaled 
eigenangles (with the same scaling by A^) of a random unitary matrix following the Haar 
measure. For the Haar measure, the law of the smallest positive point satisfies a certain 
Painleve-V differential equation (see for instance [1] for more details and references). In our 
case, one can also obtain this probability distribution as a solution of some integral equation. 
More precisely, one has the following result: 

Proposition 5.9. For all x > 0, let G{x) be the probability that the point process Too{0,jC.) 
has no point in the interval (0, x), and for all x &M., let us set: 

H{x) := l,>oX^-^G(x). 

Then H is integrable and satisfies the following equation: 

xH{x) =e I (1 - y)H{x - y) dy. (5.2) 

Jo 
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Moreover, if the Fourier transform H of H is given by 

/oo 
e-'^'^Hix) dx, 
'OO 

then it satisfies the equation: 



H{X) = H{0) exp -i9 



1 — e ^^^ — ijj, 



dyU 



(5.3) 



for all XeR. 

Proof. The probability G{x) can be expressed as follows: 



G{x) = P 



inf {Xm/Vm} > X 



m>l 



where we recall that {ym)m>i is a GEM process of parameter 9. Conditionally on yi, the 
sequence {y'^ := ym+i/i^ — yi))m>i is also a GEM process of parameter 9. Therefore, 



P 



ini {Xm/Vm} >x\xi,yi 

m>2 



P 



inf {Xm+i/l/m} > a;(i - yi) \xi,yi 

m>l 



G(x(l-yi)). 



By taking the expectation, one obtains: 

G{x)=E[l^,/y,>,G{x{l-y^))] 
= E[(l-xyi)+G(x(l-i/i))] 

-lAl/x 

9il-yf-\l-yx)G{xil-y))dy 



since the law of j/i has density 9{\ — yY ^ with respect to Lebesgue measure. Hence: 

^(^) = ^ / (^ - ^)'"'(i - ^)^(^ - ^) '^y 

where by convention, we set G{y) = for all strictly negative y. This implies the equation 
f l5.2p in Proposition 15.91 Note that for all x > 0, 

H{x) < x^-\ 

Moreover, for x > 2, < y < 1: 

H{x -y) = {x- yf-^G{x -y)<{x- yf-^G{x - 1) 

< (x - l)''-^2(^-^)+G'(a; - 1) 

<2(^-i)+iJ(a;-l), 

and by (K2\i . 

2(«-i)+6) 

H(x) < H(x-l). 

^ ^ - 2x ^ ' 

Hence, at infinity, H decreases faster than exponentially, and H, x ^-^ xH{x) are in L^. One 
deduces that H, is well-defined and differentiable. By (15. 2 p one has for all A G M: 

^'(A) = -i9H{X)K{X), 
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where K is the function defined by K{ii) = (1 — |/)lo<y<i. Therefore: 

^(A) = ^(0) exp (-19 j k{fx) dfx\ , 
which imphes 15.31 D 
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